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5 Cross-Reference to Related Application 

This application is a continuation of U.S. Serial No. 08/931,927, filed September 17, 
1997, now U.S. Pat. No. x,xxx,xxx, which was a continuation of U.S. Serial No. 08/322,313 
filed October 13, 1994, which was a continuation-in-part of U.S. application Serial No. 

j 3 08/124,098 filed by Caid et al., on September 20, 1993 (hereinafter, Caid et al), which was 

\£\ 10 abandoned in favor of file wrapper continuation Serial No. 08/561,167 now U.S. Pat. No. 

[y 5,619,709 for "System and Method of Context Vector Generation and Retrieval." All of these 

W t related applications are incorporated by reference in their entirety. 

□ 

% Background of the Invention 

i y 

: . 3 

~* 1 . Field of the Invention 

15 This invention relates to vector-based meaning-sensitive information storage and retrieval 

systems, and more particularly to an improved system and method for generating and retrieving 
context vectors that represent high-dimensional abstractions of information content. This 
invention also relates to storage, indexing, and retrieval of image data, and more particularly to a 
method and system for generating and retrieving context vectors that represent high-dimensional 

20 abstractions of information in images. 
2. Description of the Related Art 

Conventional methods of record storage and retrieval generally involve storage of all 
records word for word and then searching for key words in the records using inverted indexes. 
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The key word searches are performed by doing a complete search through all of the contents of 
the data base that contain a list of query words. Such systems have no knowledge, for example, 
that "car" and "automobile" represent nearly the same meaning, so the user must include this 
information by using a complex and difficult-to-formulate query. Some systems try to solve this 
5 problem by a built-in thesaurus, but such systems lack "meaning sensitivity" and miss many 
subtleties of meaning association, such as the fact that "car" is closer to "road" than to 
"hippopotamus". 

There is currently much research and development in the field of neural networks. A 
neural network consists of a collection of cells and connections among cells, where every 
10 connection has an associated positive or negative number, called a weight or component value. 
Each cell employs a common rule to compute a unique output, which is then passed along 
connections to other cells. The particular connections and component values determine the 
behavior of the network when some specified "input" cells are initialized to a set of values. The 



=! 15 program does in determining the behavior of a computer. 

;Jj Prior art for document retrieval includes systems using variable length lists of terms as a 

'C* representation, but without meaning sensitivity between terms. In such systems, pairs of terms 
are either synonyms or not synonyms. 

So-called "vector space methods" can capture meaning sensitivity, but they require that 
20 the closeness of every pair of terms be known. For a typical full-scale system with over 100,000 
terms, this would require about 5 billion relationships — an impractical amount of information to 
obtain and store. 

Methods have also been proposed for searching with fixed-length vectors. However, 
such methods require work on the order of at least the square of the sum of the number of 
25 documents and the number of terms. . This is impractical for a large corpus of documents or 



component values play roughly the same role in determining neural network behavior as a 



terms. 
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A document retrieval model based on neural networks and capturing some meaning 
sensitivity has been proposed. However, a search in such models requires multiplications for 
twice the product of the number of document and the number of keywords for each of a plurality 
of cycles. 

5 Koll in "WEIRD: An Approach to Concept-Based Information Retrieval," SIGIR Forum. 

vol. 13, no. 4, Spring 1979, pp. 32-50, discloses a retrieval method using vector representations 
in Euclidean space. The kernel or core used by Koll are non-overlapping documents. This 
results in rather small dimensional vectors on the order of seven values. Vectors are generated 
from the core documents based on whether or not a term appears in a document. As an 
20 alternative, Koll suggests starting with a kernel of terms which never co-occur. 
^ In addition, analysis of image subject content is a time-consuming and costly operation. 

| J J This analysis is often required for the identification of images of interest in existing image data 

^ bases and the routing and dissemination of images of interest in a real-time environment. The 
W conventional approach is to rely upon human intellectual effort to analyze the content of images. 

k 25 It would be desirable to reliably translate image data into representations that would enable a 

G 

computer to assess the relative proximity of meaning among images in a database. 

ru 

%j It would be desirable to have a computing system that can derive accurate, efficient, and 

□ manageable representations of images for later recall, retrieval, and association. 



Summary of the Invention 
20 The present invention is directed to a system and method for generating context vectors 

for use in a document storage and retrieval system. A context vector is a fixed-length series of 
component values representative of meaning or content. Geometric relationships among context 
vectors are representative of conceptual relationships among their associated items. Thus, two 
information items having similar meaning or content have similarly-oriented context vectors, 
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while items having dissimilar meaning or content have orthogonal context vectors. Similarity 
between items may be measured by calculating the dot product of the associated context vectors. 

Context vectors may be associated with words, terms, documents, document portions, 
queries, images, quantitative data, people, or any other type of information item. This use of 

5 context vectors provides a context-sensitive information retrieval, routing, and visualization 
system based on learned similarity of meaning. 

The present invention provides a context vector generation scheme that uses a neural 
network operating on a training corpus of records. Resulting vectors are relationship-based, 
formed by the proximity of words to one another in the training records. Relationships among 

10 context vectors are deterministic, so that a context vector set has one logical solution, although it 
may have a plurality of physical solutions. No human knowledge, thesaurus, synonym list, 
knowledge base, or conceptual hierarchy, is required. 

Record storage according to the present invention is performed by inputting each record, 
which may be any type of information item, in machine-readable form into a processing system. 

15 If a record is textual in nature, uninteresting words may be removed from consideration prior to 
further processing. A learning law is then applied to each word (or analogous component) of the 
record, which assigns context vectors in accordance with word proximity. The learning law 
employs a technique of "windowed co-occurrence'* wherein a fixed-size moving window is 
applied throughout the document, and words within the window (neighbor words) may exert 

20 "influence" on neighbor words in accordance with mutual co-importance. Such "influence" is 
constrained to avoid convergence, or collapse, of context vectors. Once context vectors are 
established for all components of a record, they are combined to form a summary vector for the 
record. This summary vector is then normalized. It represents the overall meaning or content of 
the record. 

25 Summary vectors of records may be clustered to reduce searching time, by forming a tree 

of clustered nodes. A centroid is computed for each node based on the average of the normalized 
summary vectors within a defined group. Thus, records with similar information content occupy 

- 4 - 5029677 



the same cluster. Information content of the cluster is represented by the centroid vector. The 
node tree may be traversed to provide a speedy method of locating relevant records in the data 
base. 

Once the context vectors are determined, records may be retrieved using a query interface 
5 that allows a user to specify content terms, Boolean terms, and/or document feedback. The basic 
searching method involves converting the user's query into a context vector (a query vector). 
Queries may also be based on example documents, selected paragraphs, sentences, or words. 
The query vector is then combined with the normalized summary vectors (or with node 
centroids) to locate the records having the closest vectors. Retrieved records may be displayed in 
20 order of vector proximity, which corresponds to relative relevance to the query. Rank ordering 
! ~ by proximity prevents information overload to the user (unlike conventional Boolean search 
! methods, where search results may include, for example, 500 documents, without any indication 

"''=4 

i-5 s 

\}{ of which documents are likely to be the most relevant). In addition, the most relevant portions of 

ii I 

•JA retrieved records may be highlighted if desired. 

"'■ci 

» 25 The system may also employ relevance feedback, whereby the user specifies which of the 

: -3 retrieved documents are most helpful. A new search may then be performed using the summary 
\| vector for the specified documents. This technique reduces the time required for searches and 
□ improves system effectiveness. 

Furthermore, the techniques of the present invention facilitate visualization of textual 
20 information by translating context vectors into visual and graphical representations. Thus, a user 
can explore visual representations of meaning, and can apply human visual pattern recognition 
skills to document searches. 

Finally, the present invention provides a language-independent information representation 
scheme. Thus, it may be used to perform English queries on foreign text for retrieval (and vice 
25 versa) without the need for prior translation or interpretation. 

The present invention is also directed to a method and system for generating context 
vectors associated with images in an image storage and retrieval database system. In this aspect 
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of the present invention, context vectors are associated with information elements, or features, 
derived by performing wavelet transformations at a plurality of points on each electronically 
stored image in the database. The transformations provide orientation-sensitive spatial 
frequencies on the images at a variety of orientation/ frequency combinations. These features are 
combined to form image feature vectors or "image vocabulary" elements analogous to words in 
text. 

A prototypical subset of feature vectors, or atoms (also called information elements), are 
derived from the set of feature vectors to form an "atomic vocabulary." In one embodiment, the 
prototypical feature vectors are derived by using a vector quantization method (e.g., self 
organization) in which a vector quantization network is also generated. 

The atomic vocabulary is used to define images in the database or any new image in 
electronic computer-readable form. As above, a wavelet transformation is performed at a 
plurality of sample points on the image to generate feature vectors representing the image. The 
generated feature vectors are mapped to the closest atoms in the atomic vocabulary using the 
vector quantization network. Thus, new images are defined in terms of the established atomic 
vocabulary. 

In one embodiment, a "stop list" of high-frequency, low-information, feature vectors is 
also generated. The stop list can be used to remove high-frequency, low-information, feature 
vectors when using the atomic vocabulary to represent images. 

In order to quantify conceptual relationships among atoms in the atomic vocabulary (and 
the images they variously represent), context vectors are employed. A context vector is 
associated with each atom in the atomic vocabulary. A learning law is applied to modify the 
context vectors as a function of the proximity of the atom to other atoms in the image and the 
frequency of occurrence of the atom in the image database. 

Once the context vectors are established, the context vectors associated with the atoms 
that define an image are combined to form a summary vector for the image. The summary vector 
represents the overall meaning or content of the image. 
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In one embodiment, summary vectors of images are stored in clusters to reduce searching 
time. Images with similar information content occupy the same cluster. In one embodiment, 
textual index terms are associated with images in the database, and are automatically assigned to 
new images. Thus, textual queries can be used to retrieve images. 

Images are retrieved using any of a number of query methods (e.g., images, image 
portions, vocabulary atoms, index terms). The query is converted into a query context vector. A 
dot product calculation is performed between the query vector and the summary vectors to locate 
the images having the closest vectors. Retrieved images are displayed in order of vector 
proximity, which corresponds to relative relevance to the query. In one embodiment, retrieved 
images are broken into sub-portions and the most relevant portions matching the query vector are 
highlighted in the image. 

Brief Description of the Drawings 

FIG. 1 A is a block diagram of the hardware elements of the present invention. 
FIG. IB is a block diagram of the training system of the present invention. 
FIG.s 2 A through 2F show an example of window definition. 
FIG. 3 is a flowchart of the training process of the present invention. 
FIG. 4 is a block diagram of the process of determining a summary vector for a 
document. 

FIG. 5 shows sample icon attributes. 

FIG. 6 is a sample display for document visualization. 

FIG. 7 shows examples of icons representing associations. 

FIG. 8 shows an example of user navigation through vector space. 

FIG. 9A is a sample display for document visualization. 

FIG. 9B is a sample display showing icon qualification. 

FIG. 9C is a sample display showing an open icon. 



5029677 



FIG. 10 is a flowchart of the retrieval method of the present invention. 
FIG. 1 1 is a flowchart of the method of document visualization. 

FIG. 12 is a chart of ranked eigenvalues derived from a set of approximately 15,000 280- 
dimensional context vectors. 
5 FIG. 13 is a chart of ranked eigenvalues derived from a set of approximately 15,000 512- 

dimensional context vectors. 
FIG. 14 is graphical representation of the hierarchical a data object definition and 

organization in a Docuverse system. 
FIG. 15 is a graphical representation of possible visual degrees of freedom for icons used 
w in the Docuverse system. 

.ttzs. 

li f FIG. 16 is a graphical representation of the coordinate selection and information display 

;■• I used in the Docuverse system. 

'*4 





FIG. 


17 


is 


a graphical representation of icon displays used in the Docuverse system. 




FIG. 


18 


is 


a graphical representation of using a "point of view" in the Docuverse system. 


15 


FIG. 


19 


is 


an example of specifying coordinates in the Docuverse system. 




FIG. 


20 


is 


an example of qualifying and selecting data in the Docuverse system. 




FIG. 


21 


is 


an example of changing the point of view in the Docuverse system. 




FIG. 


22 


is 


an example of selecting icon words and attributes in the Docuverse System. 




FIG. 


23 


is 


an example of zooming in the Docuverse system. 


20 


FIG. 


24 


is 


an example of inspecting data in the Docuverse system. 




FIG. 


25 


is 


a data flow diagram of context vector generation in the MatchPlus system. 




FIG. 


26 


is 


a continuation of FIG. 25. 




FIG. 


27 


is 


a data flow diagram of routing query generation in the MatchPlus system. 




FIG. 


28 


is 


a data flow diagram of retrieval in the MatchPlus system. 


25. 


FIG. 


29 


is 


a simplified block diagram of a computer hardware and software 



implementation of an image retrieval system in accordance with the present invention . 



5029677 




FIG. 30 is a high level flowchart of an image retrieval system in accordance with the 
present invention. 

FIG. 31 is a flowchart of one embodiment of the method for generating an atomic 
vocabulary. 

5 FIG. 32 is a flowchart of one embodiment of the method for generating a "stop list." 

FIG. 33 is a flowchart of one embodiment of the method for representing a new image 
using the atomic vocabulary. 

is a flowchart of the method for initializing context vectors, 
is a flowchart of one embodiment of the boot-strapping procedure, 
is a pictorial representation of FIG. 35. 

is a flowchart of one embodiment of the method for generating a summary vector 

is a high level flowchart of the clustering method of the present invention, 
is a flowchart of one embodiment of the clustering method, 
is a high-level flowchart of one embodiment of the indexing method, 
is a flowchart of one embodiment of the retrieval method, 
is a flowchart of one embodiment of the highlighting method, 
is a pictorial representation of highlighting as used in the present invention. 



Q 



FIG. 34 
FIG. 35 

10 FIG. 36 

FIG. 37 
for an image. 
FIG. 38 
FIG. 39 

25 FIG. 40 



■5 FIG. 41 

ru 



FIG. 42 
FIG. 43 
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Description of the Preferred Embodiments 

Section I. Context Generation and Retrieval 

Referring now to FIG. 1 A, there is shown a block diagram of a typical implementation of 
a system 100 in accordance with the present invention. The user supplies queries to system 100 
5 via input device 111. Central processing unit (CPU) 107 runs software program instructions, 
stored in program storage 112, which direct CPU 107 to perform the various functions of system 
100. In the embodiment illustrated herein, the software program is written in the C programming 
language and runs under the UNIX operating system. Each of these languages may be run on a 
variety of conventional hardware platforms. Data storage 109 contains a corpus of documents, as 
20 well as data describing context vectors. Alternatively, the corpus of documents (or other 
~3 information items) may be remotely located, with electronic links to system 100. In accordance 
=fi with the software program instructions, CPU 107 accepts input from input device 111, accesses 
*h data storage 109, and uses RAM 108 in a conventional manner as a workspace. CPU 107, data 

"i s 

^ storage 109, and program storage 112 operate together to provide a mechanism for generating 
n 15 context vectors and for retrieving information in response to queries. 

In the embodiment illustrated herein, CPU 107 can be a mainframe computer or a 
powerful personal computer; RAM 108 and data storage 109 are conventional RAM, ROM and 
□ disk storage devices for the CPU; and output device 1 10 is a conventional means for either 

I ""a 

printing retrieved information items, displaying the information on a video screen using a 
20 window-based interface system, or sending it to a database for later access. 

The preferred embodiment of the present invention is capable of handling a very large 
corpus, containing over 10 million documents. The architecture supports operation in a 
distributed data /distributed processor environment, if desired. It may be implemented on any 
platform, operating system, and user interface of sufficient power and flexibility, such as: 
25 XWindows/MOTIF; Sun/OS SunView; Microsoft Windows, VAX/VMS, and the like. 

The present invention is based upon a representation scheme using context vectors. A 
context vector is a fixed length vector having a plurality of component values that are determined 
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based on relationships between meanings of information items. Such information items may be 
words, paragraphs, queries, documents, images, and the like. In the following discussion, for 
illustrative purposes, context vectors are described with reference to words and documents, 
although many other types of information items may be similarly represented. In the preferred 
embodiment, each context vector has 200 or more component values. 

The context vector generation scheme of the present invention is designed to produce 
vectors that represent the relative proximity of meaning or content among words or documents in 
a quantitative, geometric manner. Thus, information items having similar meanings have closely 
aligned vectors, while information items having dissimilar meanings have orthogonal vectors. 
This representation scheme allows proximity of meaning to be assessed by performing a simple 
dot product (inner product) operation on associated context vectors; the higher the dot product 
result, the more similar the meanings. 

Accordingly, the absolute orientation of a particular vector in the vector-space is 
irrelevant, as long as the relative orientation (with respect to other vectors) is representative of 
relative proximity of meaning and content. In other words, the problem of finding a set of 
vectors defining relative meaning has an infinite number of physical solutions in vector-space 
(absolute orientations) but only one logical solution (relative orientations). The context vector 
generation scheme of the present invention is designed to arrive at the logical solution through a 
deterministic training method, without regard to absolute orientation of the vectors themselves. 
Context Vector Training 

Context vectors are developed for individual words or terms based on proximity to other 
words. This learning technique is performed on a training set of documents. Referring now to 
FIG. IB, there is shown a block diagram of the training system. A training text 101, stop list 
103, and phrase list 104 are provided to a preprocessor 102. Training text 101 includes a set of 
documents for training. Stop list 103 includes a list of words that are deemed uninteresting and 
are not to be considered in training (e.g., prepositions and common words). Phrase list 104 
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includes a list of multiple-word phrases that are to be treated as a single word for training 
purposes (e.g., "world series", "golden parachute", "best man"). 

Referring now also to FIG. 3, there is shown a flowchart of the training process. The 
system starts by preprocessing the documents in the training set. Preprocessing consist of several 
5 steps, including: 1) removing stop-listed words from the set of training words; 2) consulting 
phrase list 104 to locate and mark multiple-word phrases that are to be treated as a single word; 
and 3) reducing words to "stems" in order to increase the effectiveness of the training process — 
thus, "investments", "investor", and "investing" share the stem "invest" and may be treated alike. 
The set of word stems generated by preprocessor 102 is fed to learning system 105 which 
30 generates a set of stem context vectors 106 according to the method shown in FIG. 3. Each 
context vector consists of a fixed number of components (200 or more in the preferred 
embodiment). 

Learning system 105 generates stem context vectors as follows. First, initial conditions 
are assigned 303. In the preferred embodiment, initial conditions are assigned by generating a 
* 15 random context vector for each stem, consisting of components selected by zero-mean, unit- 

hU variance Gaussian random number generation. Since the system uses dot products as the 

i y 

\i measure of relationship strength, mutual orthogonality is a desirable initial condition. This is due 
□ to the fact that near-orthogonal vectors will have dot products close to zero. This near-zero dot 
product corresponds to a weak initial relationship. Assigning a random context vector provides 
20 an initial condition that approximates mutual orthogonality. As will be recognized by those 
skilled in the art, other techniques of assigning initial conditions can be employed. 

The system then starts with the first document 304 and proceeds through every document 
in the training corpus. For each document, it starts at the first word stem 305 and passes through 
the document, targeting each word stem, one by one. As each stem is targeted, the system 
25 applies 306 a learning law to the target. In the preferred embodiment, step 306 involves the 
following substeps. First, a window is defined, consisting of a fixed number of word stems 
appearing on either side of the target stem. In the preferred embodiment, the window includes 

- 12 - 5029677 




three stems on each side of the target stem, although the window can be of any size. The stems 
within the defined window are called neighbors. 

Referring now also to FIGs. 2A through 2F, there is shown an example of the window 
definition for the first few word stems 201 of a sample document. In FIG. 2 A, the target stem 
202 is "Federal" and the neighbor stems 203 are "antitrust", "law", and "simple". Window 204 
only includes three neighbor stems 203 because there are no neighbors to the left of target 202 
(since target 202 is the first word stem of the document). Next, as shown in FIG. 2B, the target 
stem 202 is "antitrust", and the neighbor stems 203 are "Federal", "law", "simple", and 
"commercial". FIG.s 2C, 2D, and 2E proceed similarly. FIG. 2F shows the more general case 
where window 204 includes three neighbors 203 on each side of target 202. (Note that "Sherman 
Act" is treated as a single word stem; this is defined in phrase list 104, discussed above). 

For each target 202, context vectors of neighbors 203 are used to "influence" the context 
vector of target 202. The relative influence of each neighbor is weighted by two factors: 1) a 
function dependent on the neighbor's position in the window relative to the target, and 2) a 
frequency function determined by the number of documents containing the neighbor stem 
(frequency). The closer the neighbor, and the lower the frequency function, the more "influence" 
the neighbor has. These weighting techniques are described below. 

The target vector Tj of a target wordj is updated using some sort of learning law. Several 
different learning laws have been developed, any of which may be used. In some learning laws, a 
weighted average of the neighbor context vectors is determined and then applied to the target 
vector Tj; in other learning laws, the influence of neighbors is determined and applied 
individually. The most effective learning law has been found to be the following: An error Ejj 
can be defined for each neighbor, representing the difference between the neighbor vector and the 
target vector. Thus: 



Eij^Nij-Tj 



(Eq. 1) 



where: 



Njj is the context vector for neighbor of the target word i of target stem j; and 
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Tj represents the context vector for target j. 

A correction value Q for target j is determined using the formula: 

ws 

q= Z -aij)Eij 

1 (Eq. 2) 

where: 

5 WS is the window size in stems; and 

q is a proximity constraint that limits the correcting effect of each error vector Eij on the 
vector Tj of targetj according to a relationship between neighbor and target,*. 

If desired, a proximity weight for each neighbor may be applied in generating the 
correction value Cj. The proximity weight causes neighbor stems that are closer to the target to 
w have a greater effect than those that are farther from the target. 

The correction value is then applied to the target vector as follows: 

3 J n r=l ' 

(Eq. 3) 



; b O where: 

'y 7 is a predefined step size, or learning rate; 

Q 

s~ 15 Fj is the total number of occurrences of stem j in the corpus; and 

M is the mean context vector for all unique stem vectors in the corpus. 
Although the above-described learning law has been found to be highly effective, 
alternative learning laws are herein presented for illustrative purposes. As will be apparent to 
those skilled in the art, other learning laws may also be applied. 
20 One alternative operates as follows. First a weighted sum Wj of the neighbor vectors for 

targetj is determined using the formula: 

ws GO) 

(Eq. 4) 
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where: 

G(i) is a Gaussian weight for the neighbor i; and 

Dj is the number of documents that contain target stem j. 

Thus, the weighted sum is inversely proportional to the number of documents that contain 

5 target stem j. The weighted sum Wj is then added to the target vector as follows: 
T NEV =T OLD +Wj 

(Eq. 5) 

Another alternative learning law, known as a "batched normalized error" law, operates as 
follows: First the weighted sum calculated in Eq. 4 is used to determine a weighted error: 

E r W r T s 

O (Eq.6) 

lift w A correction value Cj for target j is determined using the formula: 



y 



1 t=l J 



(Eq. 7) 



Then the correction value is applied to the target vector as follows: 

(Eq. 8) 

In Eq. 8, the weighted error is not normalized. Alternatively, it may be normalized, 
15 yielding the equation: 

T NElV =T OU>^£ rM 

(Eq. 9) 

One problem with these learning laws is that they may lead to collapse or convergence of 
the context vectors (overtraining) when multiple passes are attempted. The solution is to 
introduce constraints that prevent vectors from getting too close to one another. This was done 
20 in Eq. 2. The effect of a constraint is as follows. When the magnitude of the error vector is 



-15- 



5029677 



10 



ess" 

i U 



15 



20 



greater than a, the target is moved towards its neighbors. When the magnitude of the error vector 
is less than a, the target is moved away from its neighbors. Thus, convergence and collapse of 
the context vector set are avoided. 

For effective training, a should be selected in accordance with the relative "importance" 
of the target and its neighbor. If two stems are "important" and they co-occur, then the learned 
relationship between them should be strong (i.e., the dot product of the two vectors should be 
large). Therefore a should be small, permitting the two vectors to get closer to each other. In 
order to vary a according to the importance of both the neighbor and the target, the following 
technique is employed. 

Importance is determined by the frequency of occurrence of the stem in the corpus. The 
importance metric Ij for stem j varies from a predefined lower bound B to 1.0: 



Importance is determined by the equation: 



(Eq. 10) 



/ 



/y=B + (l-fi) 



log 



1- 



{nDj) ' 



\ 



( TND ) 



log 



/ 



(Eq.ll) 



where: 

NDj is the number of documents that contain stem j; and 
TND is the total number of documents in the corpus. 

From this equation, it can be seen that a stem appearing in every document will have an 
importance equal to B, while a stem appearing in only one document will have an importance 
equal to 1 .0. The greater the frequency of occurrence, the less important the stem. 

The value of a for a particular pair of stems i and j is determined by the equation: 
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Thus, the greater the co-importance of the two vectors, the smaller the value of a. The 
value of Of is bounded as follows: 

(Eq. 13) 

It can be seen, therefore, that the value of Of determines how close any neighbor vector 
can get to the target vector. The value of Q! determines the minimum angle between the two 
vectors, and thus, the maximum dot product between them. 

Referring again to FIG. 3, whichever learning law is used, the system then checks 307 to 
see whether there are more stems to be targeted in the current document. If so, it targets the next 
stem 308 and returns to step 306. 

Once the context vectors for all stems of the document have been targeted, the system 
determines 329 a summary vector for the document as a whole. This vector is representative of 
the overall content or meaning of the document. It may be generated by simply adding the 
context vectors of all the stems in the document and normalizing the result. Alternatively, stems 
may be weighted according to their frequency; the context vectors of stems that occur less 
frequently in the corpus as a whole are weighted more heavily when the summary vector for a 
document is calculated. Other techniques of determining a summary vector from a set of stem 
context vectors may also be employed. Referring now also to FIG. 4, there is shown a block 
diagram of the process of determining a summary vector for a document according to the 
preferred embodiment. Context vectors 106 are weighted according to an inverse frequency 
weight 401 and combined with a vector summing operation 402. The result is normalized 403 to 
produce a normalized summary vector 404. 

Referring again to FIG. 3, if there are more documents to be processed, the system goes to 
the next document 311 and repeats steps 305 through 309. 

Target vectors are not actually updated until the end of an iteration (one pass through the 
corpus). This prevents flip-flopping and unwanted feedback effects as vectors move towards 
their neighbors. All corrections are summed during the iteration and applied at iteration 
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boundaries, in a vector update 312. After all targets Tj are updated, the mean M is also updated 
313. 

The system then checks 314 whether additional iterations are required for the corpus. 
The number of iterations to be performed for a document depends upon some kind of predefined 
stopping criterion that may incorporate speed, stability, or other concerns. If additional iterations 
are required, the system returns to step 304. 

As will be apparent to those skilled in the art, many variations on the above training 
techniques are possible, particularly where the training data are non-textual in nature (such as 
imagery, sound, video, and the like). For illustrative purposes, a more generalized description of 
a context vector generation technique is given in Development of Context Vectors by Singular 
Value Decomposition , below. 
Context Vector Storage and Indexing 

Once the summary vectors have been determined, they are stored. Storage of the 
normalized summary vectors can be arranged to further reduce searching time by creating cluster 
trees. An initial parent node at the top of the tree indexed as level 0, node 1, initially contains all 
of the normalized summary vectors in the data base. A series of child nodes each branching from 
the initial parent node is created at a next level of the cluster tree. A centroid-consistent 
clustering algorithm is used to distribute the summary vectors among the series of child nodes. A 
group of clusters is centroid-consistent if every member of every cluster belongs to the cluster in 
the group with the closest centroid. A centroid is determined by taking, for each dimension, an 
average of the component values from all of the context vectors within the group. One popular 
centroid-consistent clustering algorithm is convergent k-means clustering. Convergent k-means 
clustering can be performed as follows: 

1 . Begin with any initial partition that groups the vectors into k clusters. For 

example, take the first k summary vectors as single element clusters and set the initial 

value of the centroid of each cluster to equal its member vector. Assign each of the 
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remaining summary vectors to the cluster having the nearest centroid. After each 
assignment, recompute the centroid for the cluster which gains a vector; 

2. Take each summary vector in sequence and compute its distance from the 
centroid of each of the k-clusters. If the vector is not currently in the cluster with the 

5 closest centroid, move the vector to that cluster and update the centroids of the clusters 

that gain or lose a summary vector; 

3. Repeat step 2 until convergence is achieved, that is until a pass through all of 
the summary vectors causes no new assignments. 

Since convergence may be rather time-consuming to achieve, the clustering algorithm can 
10 be simplified by limiting the number of repetitions of step 2. After a given number of repetitions, 
"f such as 99, the centroids can be frozen. Then, one or more passes can be made through all of the 
\p summary vectors, distributing the vectors to appropriate clusters, but without updating the 

centroids. While, using this approximation, the centroids will no longer be exact centroids, the 
approximate centroids will be sufficient for the use of the present invention. It is not necessary 
« 15 to the present invention that the centroids be precise; however, it is preferable that the clusters be 
us centroid-consistent. The last pass through the summary vectors guarantees that the clusters are 

%j centroid-consistent with the approximate centroids. From herein, "centroids" as used in this 

Q 

r 3 application shall mean approximate centroids — in other words, a centroid sufficient to establish 
centroid-consistent clusters. Each node is identified by its centroid for use in the searching 
20 process. 

Alternatively, other centroid-consistent clustering algorithms may be employed. Such 
alternatives are described, for example, in "Self-Organized Formation of Topologically Correct 
Feature Map"; T. Kohonen, Biological Cybernetics, vol. 43, p. 50-69, 1982; and "Performance 
Evaluation of Self-Organized Map Based on Neural Equalizers in Dynamic Discrete - Signal 
25 Detection"; T. Kohonen, et al., in Artificial Neural Networks; and "The Self-Organizing Map", 
Tuevo Kohonen, Proceeding of the IEEE, Vol. 78, No. 9, September, 1990, all of which are 
incorporated herein by reference. 
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In forming a next level of clusters, the nodes in the level above become parent nodes to a 
set of child nodes below. Only the summary vectors assigned to a parent node are used in the 
clustering algorithm to form the child nodes which branch from that parent. This is repeated 
across the entire level of parent nodes and on subsequent levels so that fewer and fewer context 
5 vectors are assigned to the child nodes on each lower level. The nodes form a tree pattern in 
which each node branches from a node in the level above. Thus, each summary vector is 
assigned to a node on each level of the cluster tree. Each node has a centroid. The bottom-level 
node assignments for each summary vector and the centroids for each node are stored for use in 
the search and retrieval algorithms. On the bottom level of the tree, each node points to each 
10 normalized summary vector assigned to it. The nodes on the bottom level may be referred to as 
~ buckets. 

\*l Once a cluster tree has been set up, it is a simple matter to add a new record summary 

!%l vector to the tree. The initial branches of the tree are examined to find the closest centroid. The 

ru 

j»U summary vector is assigned to the node with the closest centroid. Then, the branches from that 

« 15 node are examined for the closest child node centroid, and the process is continued until a bucket 
v3 is reached. The new record is then assigned to the bucket with the closest centroid of those 

;:: ; 

: H 

y buckets branching from the node on the previous level to which the summary vector was 
q assigned. The centroids themselves are not changed. This action preserves centroid-consistency 
of the clusters. If a bucket gets too large, the summary vectors on the bucket can be divided into 
20 subclusters on a subsequent level. 
Retrieval 

Referring now to FIG. 10, retrieval of context vectors that have been stored according to 
the above-described tree technique proceeds as follows. The user makes an inquiry using a set of 
words or terms, or by specifying one or more records (or documents) for which similar records 
25 are sought. The words, terms, and records may be weighted if desired to designate which are 
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most pertinent to the results being sought. After appropriate weighting, context vectors for the 
words, terms, and records are combined by addition to produce a single query vector 1002. 

The query vector is then compared with each summary vector in the data base by 
performing a dot product vector operation 1003. Lists of these results are accumulated 1004, and 
5 the summary vectors resulting in the highest dot products are considered most relevant to the 
search. If desired, documents can be ranked 1005 in order of dot product magnitude to form a 
list indicating relevance. The ranked list is then output to the user 1006. Rank ordering by 
proximity prevents information overload to the user (unlike conventional Boolean search 
methods, where search results may include, for example, 500 documents, without any indication 
w of which documents are likely to be the most relevant). 
: f By using the above-described cluster tree storage mechanism, the searching task can be 

I; 1 j greatly accelerated. The query vector is used to identify the summary vector that is closest to the 
jj* query vector (by dot product computation). The search is performed using a depth first tree walk. 
yJ A branch is followed down the tree, taking the node at each level having the centroid closest to 
» 15 the query vector. The search proceeds down the tree until a bottom level node (bucket) without 
;J3 children is reached. Each of the summary vectors in the bucket is compared with the query 
y vector (again by dot product computation) to identify the closest summary vector. 

s 

13 Before a subsequent node in the depth first tree walk is checked for a closest vector, first 

it is determined whether the node can be completely pruned. A node is pruned if it is not 
20 possible for a closer normalized summary vector to be assigned to the node than the closest 

normalized summary vector found so far without violating centroid-consistency. Suppose we are 
examining a node with centroid C for pruning. If C is the centroid of any sibling node, then if it 
is true that any vector closer to the query vector Q than V (closest vector found so far) must be 
closer to C than C\ then we may prune the node with centroid C as well as any nodes branching 
25 therefrom. This may be computed by comparing the distance between C and C with twice the 
sum of the distance between C and Q and the distance between Q and V. If the distance between 
C and C is greater, then the node with centroid C (and its descendants) may be pruned. If not, 
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the formula is repeated for the remaining sibling nodes, since any one of them may permit 
pruning to proceed. If none of the sibling nodes achieve pruning of the node, then the search 
continues through the node with centroid C and down into the subsequent level if there is one. 
By using the pruning formula, a node can be pruned when any vector closer to the query vector 
than the closest vector V must be closer to the centroid C than to the centroid C\ Therefore, that 
vector could not be assigned to node C or else it would violate centroid-consistency. If this is a 
bottom node, then all of the summary vectors on the node must be checked to determine whether 
any are closer than the closest vector found so far. If a closer summary vector is found, it will 
then become the closest summary vector being remembered. Thus, bottom nodes are thoroughly 
searched if not pruned. The search continues in a depth first tree walk, pruning off entire 
branches when possible. These searches continue through the tree until all branches have either 
been checked or pruned. After the entire tree has been searched, the closest summary vector has 
been identified. The record associated with the summary vector can be retrieved. 

Other known techniques for node pruning, including linear programming techniques, may 
also be employed. 

The system may also employ relevance feedback, whereby the user specifies which of the 
retrieved documents are most helpful. A new search may then be performed using the summary 
vector for the specified documents. This technique reduces the time required for searches and 
improves system effectiveness. 

In addition, the most relevant portions of retrieved records may be highlighted if desired. 
This is done by dividing each retrieved records into a number of sections, representing chapters, 
paragraphs, or other components. A summary vector is generated for each section, based on the 
word stems in that section. Dot product computation of the section summary vectors with the 
query vector is then performed to isolate those sections that are most relevant to the query. The 
selected sections are then displayed using some distinguishing visual attribute (bold, larger type, 
a different font or color, an enclosing box, and the like). Thus, the user is able to quickly locate 
the portions of the document that are most relevant to the query. 
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One of the possible applications of the above-described system is in the area of 
conventional ICD9 codes that are commonly used to describe medical procedures. For example, 
context vectors could be developed to represent medical procedures and their associated ICD9 
codes. Then, when additional information is needed for a medical procedure, a query vector 
could be formulated to retrieve procedures and codes that are relevant to the current procedure. 

Another application of the information retrieval system described above is the automated 
coding of text documents according to a defined index of terms. For example, the Wall Street 
Journal uses an index of approximately 150 terms to code each article. These terms are assigned 
by human editors. The information retrieval system described above can be used to emulate the 
performance of the human editor in assigning such index terms, in the following manner: 

1 . Build context vectors for words using a sample of text. 

2. Collect a set of documents that have been indexed by human "experts" 
(e.g., editors in the case of the Wall Street Journal), called the indexed collection, and 
generate context vectors for these documents. 

3. Generate a context vector for the new document to be automatically 
indexed. 

4. Compare the context vector of the new document with the context vectors 
of all documents in the indexed collection, and identify the best matches (perhaps the 10 
best matches). 

5. Produce a list of the index terms of each of the best matches and assign a 
weight to each term that is proportional to the degree of match such that better matching 
indexed documents have larger weights than indexed documents that have do not match 
as well. 

6. For each unique index term, generate the index term score by adding the 
weights of each occurrence of that index term in each of the best matching index 
documents. 
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7. Sort the list of unique index terms according to the index term score and 
assign to the new document those index terms at the top of the list. 
See HNC MatchPlus Functional Specification section below, which provides a functional 
specification, including module and data format descriptions, for a preferred context vector 
5 generation, storage, and retrieval system according to the present invention. 

Document Visualization 

Another useful application of stored summary vectors is in the area of visualization of 
document content. Context vectors provide a mechanism by which meaning and content of 
documents can be represented in a visual form, allowing a human observer to take advantage of 

10 visually-oriented pattern recognition skills to find documents that are of interest. 

Once vectors have been established using the above-described methods, they can be 
represented visually using any of a number of techniques. The preferred embodiment provides a 
color graphics visual representation of a set of documents on a computer screen, such as a color 
graphics workstation or a PC or PS/2 computer equipped with a graphics board. It operates using 

25 software written in the C programming language and runs under the UNIX operating system. 
Essentially, the summary vectors for documents and other information items are displayed in a 
pseudo-orthogonal display having axes corresponding to query terms (the query terms need not 
actually be orthogonal to one another). In addition, the items may be displayed with visual 
attributes representing relative similarity of meaning with other query terms. 

20 Referring now to FIG. 1 1, there is shown a flowchart of the method of document 

visualization of the present invention. The method is described with regard to documents, 
although it may be applied to any type of information items. First, the system accepts 1 102 a 
query from the user. This query may be in the form of a word, term, phrase, document, or other 
form, as discussed above. Then, the query is separated 1 103 into a number of components. This 

25 separation may be performed manually by the user, or some automated means may be used, to 
generate components that maximize information content for the display. For example, query 
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components may be selected in alignment with the principal components of the document set 
covariance matrix. These are obtained by considering the data object's context vectors as points 
in feature space. These points form a cloud, with one point for each object. Principal component 
analysis fits the best-fitting ellipsoid to this cloud, based on root-mean-squared analysis. Query 
terms corresponding to the longest perpendicular principal axes of this ellipsoid are selected as 
the principal components of the data set. 

For each component, a display coordinate (axis) or visual attribute is assigned 1 104. For 
example, if five components are identified, the first three may be assigned to the X,Y, and Z 
axes, and the remaining two may be assigned to the visual attributes of color and texture. Any 
number of axes and attributes (including visual and non- visual attributes) may be identified and 
assigned to components, although an excess of axes or attributes may lead to a confusing display. 

A context vector is determined for each query component 1 105 by the method previously 
described. Then, for each document to be represented in the display, a set of dot products is 
computed 1 106 using the summary vector of the document with the context vector for each query 
component. The resulting set of dot product results for each document specifies coordinates and 
visual attributes for a representation of that document. Once such information has been 
developed for all documents to be displayed, the results are supplied 1 107 to a display engine 
capable of on-screen display of the icons. Thus, a set of icons are shown on a display screen, 
having on-screen positions and attributes corresponding to the dot product results. 

The display engine may employ known artificial reality image generation technology to 
portray each document as a three-dimensional icon with a specific shape, size, color, texture and 
movement projected into a higher dimension context vector space, in accordance with the dot 
product results previously determined. In addition, coded information about a document (such as 
the author or the date of publication) can also be represented. 

Many variations on the above-described visualization scheme are possible. If desired, the 
above technique may be employed using only display coordinates; all icons will then be 
displayed having uniform visual attributes. Alternatively, icons having visual attributes such as 
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color, size, and the like, could be displayed without using a positioning scheme. Thus, a 
relatively compact display could be generated, without a need for high-powered processors to 
generate the artificial reality display described above. In one embodiment, each icon contains 
one or more small images of thermometers, each thermometer indicating the degree of 
5 correlation (dot product result) with a particular concept. 

The display engine of the preferred embodiment is a high-level graphics software^ 
interface such as the Programmer's Hierarchical Interactive Graphics System (PHIGS). Other 
display engines may be used as well. PHIGS and other systems are described below, as well as 
in the following publications which are incorporated herein by reference: Hill, F.S., Computer 
w Graphics, Macmillan, New York, 1990; Kessener, L.R.A., Data Structures for Raster Graphics , 
: Z Springer- Verlag, Berlin, 1985; Foley, J.D., and van Dam, Fundamentals of Interactive Computer 
^ Graphics , Addison- Wesley, Reading, MA, 1983. 

^ The description of a three-dimensional scene used as input to PHIGS is simply the 

'U 

*U definition of each individual object to be displayed, expressed as a set of linked polygons located 

= 25 in a fixed 3-dimensional coordinate system, with each polygon having specified light reflectivity 

i3 properties (color, specularity, texture, etc.). The polygons make up the visible exterior surfaces 
'U 

■•■4 of the objects to be displayed. PHIGS handles the lighting of the objects and the calculation of 
=5 their appearance to the user from a particular vantage point. 

The use of such visual display techniques allows a user to view large groups of 
20 documents simultaneously in a multi-attribute space. The display of the present invention 
simultaneously shows the user all of the attributes of each data object for a large set of data 
objects. 

Referring now to FIG. 5, there is shown a set of examples of six different visual attributes 
501-506 for icons: size 501, shape 502, color 503, distortion 504, orientation 505, and motion 
25 506. Many other examples are possible, including additional visual characteristics as well as 
sound. 



-26- 



5029677 



Referring now to FIG. 6, there is shown a sample display having three axes 601, 
corresponding to the terms "money laundering," "drugs," and "terrorist." Referring now to FIG. 
9 A, there is shown a sample display with axes 601 and clusters of icons 901 positioned in the 
coordinate space defined by axes 601. 
5 Referring now to FIG. 7, there is shown an example of icon display. In this example, the 

user has specified two icon words: "Sendero Luminoso", having the attribute of size; and 
"Simon Bolivar", having the attribute of distortion. Thus, large icons as shown in box 701 
indicate a strong association with Sendero Luminoso, while small icons as shown in box 702 
indicate a weak association with Sendero Luminoso. Similarly, distorted icons as shown in box 
w 703 indicate a strong association with Simon Bolivar, while undistorted icons as shown in box 
; s :f 704 indicate a weak association with Simon Bolivar. 
[ In the present invention, each individually resolvable icon is portrayed as an easily 

identified object in its proper position, and possessing its assigned attributes. Icons are displayed 
W in simulated three-dimensional space, perspective, and hidden line removal. By means of simple 

15 mouse commands, the user is able to navigate through the three-dimensional projection of the 
i k 0 higher dimensional context vector space. A user-selected window is available to show the entire 

ru 

%j vector space, as well as the position and orientation of the user's current viewpoint. The position 
i«5 and orientation can be changed in response to user commands. These operations are performed 
using conventional computer graphics and artificial reality techniques. 
20 Referring now to FIG. 8, there is shown an example of user navigation. Starting view 

801 includes three objects 804 positioned with respect to three axes 805. The user selects two 
operations 803, a rotation and a translation, resulting in ending view 802. In ending view 802, 
axes 805 have moved, and objects 804 have been repositioned accordingly. 

When an icon is located at too great a distance from the user's position to accurately 
25 represent all of its characteristics, it will be shown as a point of light. If a large number of icons 
are located close to one another, they may be shown as a cloud or a shaded region. 
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If desired, the user may specify data object qualification parameters to help reduce visual 
clutter and information overload. One method of qualification is to allow the user to specify 
Boolean parameters, and only display icons which match the specified parameters. Alternatively, 
the display may be limited to the top-ranked documents resulting from a context vector query 
performed as described above. Referring now to FIG. 9B, there is shown an example containing 
axes 601 and one cluster 901 of icons that has been qualified. Other icons outside cluster 901 are 
not displayed. 

Another way to reduce information overload is to provide hierarchical organization of the 
icons. The user selects an icon for examination of subordinate objects, and specifies a "zoom in" 
command. When the system zooms in on an icon, all icons representing other documents are 
erased from the display. New "sub-icons" are introduced, representing sections, chapters, and/or 
paragraphs of the selected document. These sub-icons are displayed in the same manner as 
icons. 

Responsive to some command, such as a double-click, associated with one of the icons or 
sub-icons, the associated document (or portion of a document) may be displayed for perusal by 
the user. Referring now to FIG. 9C, there is shown an example of a text window 910 
superimposed on a display having axes 601 and an icon 909. Window 910 contains the text of 
the document associated with icon 909. The user can move or scroll, adjust window size, and 
close the window as desired, using conventional window manipulation techniques. 

Development of Context Vectors™ by Singular Value Decomposition 

A NEW APPROACH TO BUILDING MUTUAL SIMILARITY OF USAGE 

RELATIONSHIP REPRESENTATIONS 
1 . Introduction 

Context vectors are a new and powerful approach to the representation of mutual 
similarity of usage (and/or meaning) between objects in a large data corpus. Context vectors are 
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being used for document retrieval and routing and for visualization of the contents of groups of 
documents. 

Beyond these proven uses, it is also believed that context vectors can be used for retrieval 
and routing of foreign language databases and for the construction of content-addressable image, 
sound, and video databases. Applications of context vectors to the detection and localization of 
interactions between outlaws (e.g., international criminals, gangs, and terrorists) are also being 
considered. 

The technical approach discussed here builds upon a number of existing concepts, ideas, 
and results. We have combined and improved upon these important past works to achieve 
something qualitatively greater. However, at least two critical elements of our approach are 
entirely new— namely, the idea of context vectors and the observation that proximate co- 
occurrence is cognitively central to human mental association, and therefore applicable to the 
assessment of similarity of usage or meaning for almost all raw data— text, imagery, sound, video, 
etc.— see Section 2 for details. 

This White Paper addresses a new possible method for generating context vectors much 
more efficiently and more accurately than the current method allows. This new method employs 
principal component analysis, implemented using Singular Valued Decomposition (SVD) and 
other related techniques. 

To provide background for the discussion to follow, the context vector concept is defined 
and explained in the next section. Following this, Sections 3 through 6 discuss the proposed SVD 
core context vector generation method and its application to text retrieval and routing. Section 7 
discusses the use of context vectors in building a multi-language text retrieval and routing 
system. 

2. Context Vectors 

In this section the concept of a set of context vectors is defined. This definition is general 
and can be applied to text, imagery, sound, and video. 
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The fundamental idea behind the context vector concept is the observation that many 
large bodies of real-world data have the property that the statistics of proximal co-occurrence of 
their basic elements are fixed. Further, frequent proximal co-occurrence of two basic elements 
can, in almost all cases, be interpreted as implying a strong similarity of usage or "meaning" 
5 between these elements. For example, in all human languages, in both spoken and written form 
(here basic elements are words), elements that frequently appear near one another have a strong 
associational linkage that could be called "similarity of usage or meaning". Similarly, still or 
video image elements that frequently appear near one another (with "near" taking on a 
spatiotemporal meaning in the case of video) can also be taken to have a similarity of usage or 
10 meaning. The utility of these observations for data retrieval and routing probably stems from the 
fact that the associations formed by the human mind are themselves derived, at least in large 
measure, by training processes driven by common proximal co-occurrence (although, of course, 
many other things are learned about such data as well). In effect, proximal co-occurrence can be 
viewed as a first-order approximation to the associational processing of the human mind. What 
25 has been discovered during the TIPSTER Program work is that this first-order approximation 
! *D may be more than sufficient for carrying out high-performance retrieval and routing. 

ru 

%i Context vectors are the quantifications of the statistics of proximal co-occurrence. We 

□ 

H now define them. 

We assume the existence of a large data set (text, imagery, sound, video, etc.). By large it 
20 is meant that if the statistics of the structures we will study are determined using a (randomly 
chosen) sizable subset of the database that these statistics will not, with very high probability, 
change significantly if re-measured on the whole data set. 

We assume that our data set is densely endowed with what we will call elements, each of 
which belongs to one of N classes. For example, in text, the elements are words that belong to a 
25 designated lexicon (other words are ignored). Each lexicon word in the database belongs to one 
of N classes (stems). In imagery, the elements might be the objects designated by an automated 
attentional focusing system. These elements would each be assigned to a single class (where the 
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classes are self-defined as distinct "clusters" of range/azimuth/elevation/background-insensitive 
feature vectors). 

For convenience, we will assume that each element A in the database is numbered with a 
unique integer index i. We will refer to the i element of the database as A,-. The number of the 

5 class to which A,- belongs will be denoted by c ( -, where 1 ^ ^J. The set of indices of all elements 
belonging to class K, 1 ^1, will be denoted by Sg. 

Each time an element appears in the database other elements typically appear "near" it. 
For example, when the word "context" appears in this document the word "vector" often appears 
next to it. It is these proximate co-occurrences that we exploit. To carry out this exploitation we 

w must be given co-occurrence examples. Each such co-occurrence example consists of a base 
element A, along with a set of other elements found in close proximity to A,-. For each such 
proximal element Ay in a co-occurrence example with base element A/ a proximity weighting d,y 
is defined. The proximity weighting lies between 0 and 1, with greater weights representing 
greater proximity or nearness. Note that each database element can appear as the base element in 

15 at most one co-occurrence example. Some elements may not appear as base element any co- 
occurrence examples because we might not have all possible co-occurrence examples available 
to us. Finally, if we are given a co-occurrence example with base element A/ in which Ay appears 
with proximity weighting d,y, then we assume that we will also be given a co-occurrence example 
with base element Ay in which A, appears with the same proximity weighting dy, =d,y. In other 

20 words, we assume that proximity weightings are symmetric. 

Given a large set of co-occurrence examples, we then form the square mutual co- 
occurrence matrix R=(r MV ), where: 



with the averages being taken over all co-occurrence examples available for training the system. 
25 Note that since we have assumed that for each d,y we have an equal dy/ in the training examples, R 
is automatically symmetric. Also note that since d,i can be assumed to be 1 (each element co- 



r MV =AVERAGEf^J4 } ^0 and i eS u and j eS v } 9 
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occurs with itself with proximity weighting 1), that R has its diagonal elements all equal to 1. To 
support the work that follows we note that because R will typically be very sparse and because 
almost all of the off-diagonal elements will be less than 1, it is not unreasonable to make the 
following assumption: 
5 Assumption 1 : The matrix R is assumed to be positive definite. 

As shown by Strang [20], this is equivalent to saying that the determinants of the upper- 
left, diagonal-centered square submatrices of R all have positive determinants. This can normally 
be expected to be true in the cases of text, sound, imagery, and video (and can be tested for any 
specific R). 

10 Given the above preliminaries, we are now ready to define context vectors. 

; s ~ Definition 1 : A set of vectors 

J {w lt w 2 , . . . , w N }dR n , 

TU where n is an integer with n»l, is called a set of context vectors for a data set mutual co- 

us 

J occurrence matrix R iff 

□ 15 W, »Wy =T,y 

ry for all i, j=l, 2, . . . , N. 

!«3 The meaning of this definition is that the context vectors (which are all unit-length 

: ^ vectors, since w, »w/ =r,i =1) have directions that represent the mutual co-occurrences of the 
element classes. Note that since the w, are unit vectors that w, # w 7 =r,y =cos(.,y), where is the 
20 angle between w, and w 7 . As mentioned above, the context vectors encode the relative usage or, 
in some crude sense, the relative meaning of the element classes. Note that if we were to rigidly 
rotate the entire context vector set in R n that we would get an equivalent set of context vectors 
(since it is only the mutual angles between the vectors that matter). Thus, context vectors are not 
unique. 

25 Now that we have precisely defined what a context vector is, some questions arise. First 

and foremost, given a mutual co-occurrence matrix R, do a set of context vectors exist? Another 
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question is: given R, can we create a set of context vectors from it? Finally, how small can n be? 
These are some of the questions we examine in the next two sections. 
3. Context Vectors Exist And Can Be Constructed 

Given a mutual co-occurrence matrix R it is easy to show that a set of context vectors 
5 must exist for it. 

Theorem 1 : Given a symmetric, positive definite mutual co-occurrence matrix R, there 
exists a set of context vectors in R^ for it. 

Proof: If we define the NXN matrix W to be 

W=[W/, W 2 , . . . , Wyv ], 

10 where the NX1 w^ vectors are the columns of W, then the condition that the context vectors must 
I*- meet (from Definition 1) can be re-expressed as 

!;[] W r W-R. 

l J= Thus, we must show that there exists such a matrix W. To do this all we need note is that 

* y 

w since R is symmetric and positive definite we can use Gauss decomposition to reexpress this 
«_ 15 matrix as 

j s 

<S R=QDQ r , 
IV 

\j where the columns of the orthogonal matrix Q are the unit-length eigenvectors of R and the 

.:ss. 

p matrix D is diagonal, with the eigenvalues of R as its diagonal entries (see Strang [20] for 

details). If we then take the square root matrix of D (namely, the matrix D 1/2 which has each of 
20 its diagonal elements equal to the square root of the corresponding diagonal element of D) and 
rearrange terms we get 

R=QDQ r =QD //2 D //2 Q r =(D I/2 Q T ) T (D //2 Q r ). 

Thus, we can take W to be D y/2 Q r . Thus, for any symmetric, positive definite R there 
exists a W in 

25 The upshot of Theorem 1 is that if we are willing to let n=N (which, for most 

applications, will be a very high dimensional space), then we can always find a set of context 
vectors in R* no matter what R is. Namely, all we need to do is calculate the eigenvalues and 
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eigenvectors of R and then use them to construct W=D Q . While this indeed gives us a set of 
context vectors, they are not really ones that we want. Constructing desirable sets of context 
vectors in R* where n«N is the subject of the next section. 
4. SVD Construction of Practical Context Vectors 
5 While the formula W=D //2 Q r for constructing context vectors will work, it essentially 

leaves us stuck with using N dimensions for representing these vectors. For a practical 
application this would typically be a debilitating disadvantage in terms of computer memory and 
computational burden. For example, in a text data access system we might have N=250,000 (the 
number of stems in the corpus). In this section we demonstrate a method for producing context 
10 vectors in a (typically) much lower dimensional space that can perform approximately as well as 
-3 "full-dimensional" context vectors. In fact, we provide a formula for calculating the error induced 
j«P by reducing the context vectors to n dimensions from N dimensions. 

fU The basic idea is to employ the Singular Value Decomposition (S VD-see Strang [20] for 

i y 

Ly details) in which any arbitrary MXN matrix S can be written as 
/ 15 S=PD //2 Q r , 

;.n where P is an MXM orthogonal matrix (with columns given by the unit-length eigenvectors of 
vf S S), Q is an NXN orthogonal matrix (with columns given by the unit-length eigenvectors of 
}t SS 7 ), and D //2 is an MXN matrix with the square roots of the r non-zero eigenvalues of S r S (and 
SS r !~they have the same non-zero eigenvalues, all of which are positive and real) as the entries 
20 in its upper left "diagonal" and with zeros in all other positions. The orders of the eigenvectors in 
P and Q must match the order of the corresponding eigenvalues in D. Note that the SVD is not 
unique (for example, we can rearrange the orders of the eigenvalues and eigenvectors in D, P, 
and Q). From now on we will arrange D so that the eigenvalues are listed from largest to 
smallest, starting with the largest eigenvalue at the upper left corner of the matrix and moving 
25 down the diagonal to the right 

In our case we will be interested in applying the SVD to yield a lower-dimensional set of 
context vectors from our matrix R. To do this, we first note that the SVD expansion depends only 
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on the properties of the matrices S r S and SS T . Since we want W r W=R, and since R is symmetric, 
we can identify R with both S r S and SS T . Thus, the eigenvectors of S r S and SS r will be the 
same. So, in this SVD case, we will have P=Q. Thus, we can write W as 
W=QD //2 Q r , 

5 with Q and D taking on the same meaning as in the previous section. It may seem odd that this 
construction of a W is so close to, and yet clearly different from, that provided by Theorem 1 . 
However, as we noted above, W is by no means unique. To see that this SVD-derived W is 
indeed acceptable, note that 

W r W=(QD y/2 Q r ) r QD I/2 Q T =QD y/2 Q r QD //2 Q r =QDQ r =R, 
w where we have exploited the property that Q r Q=QQ r =I (the identity) for any orthogonal matrix. 

The value of using W=QD //2 Q T instead of W=D //2 Q r becomes clear when we expand 
I this new W in terms of the columns of Q. To make the notation clearer, we let 
[H Q=K] = [u/, u 2 , . . . , u N ], 

where the NX1 vectors u^ (the unit-length eigenvectors of R) are the columns of Q. Using 



-4 



25 W=QD //2 Q T we can then write 



S N in N 

Wi= Z (kj u Uij) uj= £ 

j=l 7=1 



a 



where the (necessarily non-negative and real) eigenvalues .y, .2, . . . are listed in descending order. 
By virtue of this formula, we can see that each w, vector is expressed as a weighted sum of a 
fixed set of N vectors (the u 7 ). Because of the nature of the matrix R 5 it is reasonable to expect 

20 that many of its eigenvalues j will be close to zero. This is true of almost all real-world matrices 
constructed in a like manner, and we certainly expect it to be true in this case (see below for 
some experimental results supporting this supposition). Because of this, and because the up- 
values are all small (recall that Q is orthogonal, and thus the u,y are the components of unit 
vectors), the .,y scalar coefficients multiplying the u y vectors in the w, sum become progressively 

25 smaller as j increases (regardless of which particular w, vector we are considering). Thus, it may 

- 35 - 5029677 



be reasonable to truncate these sums after r terms, where Kr<N. In this case, we get the 
approximation 

r 

Wi* Z ctijuj. 

Note that the u y - vectors are an orthonormal basis for R^. However, we have now 
discarded all but r of them. So this is an orthogonal basis for R r . Thus, we can think of our w, 
vectors as belonging to an r-dimensional Euclidean space. If we reexpress our w,- vectors in terms 
of the uy basis vectors (i.e., taking the first coordinate of this new representation to be the u 7 
component, the second to be the u^ component, and so on) we get 

wr=((a a , a^, . . . , oLi r ) T 3Vf. 

The error introduced by using these r-dimensional approximations v, can be quantified. 
Specifically, the errors in the inner products will be given by 

N 

k=r+l 

In conclusion, we have shown that a set of lower-dimensional context vectors can be 
constructed by applying the SVD to the mutual co-occurrence matrix R. The error introduced 
into the context vectors by reducing the dimension from n=N to n=r is quantifiable and can be 
controlled as required by the specific problem. 

5. Applicability of the SVD Method 

The first question that naturally arises in connection with the SVD method proposed 
above is how rapidly the eigenvalues of the co-occurrence matrix R actually fall off in a real- 
world situation. This will determine both the practical dimensionality of the context vectors and 
the utility of the SVD method. This question is at least partially answered by FIGS. 12 and 13. 
These figures were created by calculating the eigenvalues of the matrix 
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T=WW r , 

where, in this case, the W matrix is defined by a set of context vectors obtained using our current 
random-initialization and bootstrapping context vector generation method (the columns of W are 

T T 

the context vectors). As is well known, the matrices W W and W W have the same non-zero 
5 eigenvalues. However, W W is an NXN matrix (where N is the number of stems—approximately 
15,000 in these examples), whereas W W is an nxn matrix, where n is the dimension of the 
context vectors (n equals 280 and 512 in the examples of FIGS. 12 and 13, respectively). 
Therefore, we calculated T=W W r . 

FIG. 12 Ranked eigenvalues (largest on the left to smallest on the right) for T=WW r 
w derived from a set of approximately 15,000 280-dimensional context vectors (i.e., N 5,000 and 
! ~" n=280). Note that the smallest eigenvalues are not much smaller than the largest eigenvalues. 

; This suggests that the dimensionality of this context vector space cannot be reduced further 
| E H (perhaps increasing the dimensionality would improve somewhat). 

W FIG. 13 Ranked eigenvalues (largest on the left to smallest on the right) for T=W W r 

s 25 derived from a set of approximately 15,000 512-dimensional context vectors (i.e., N=4 5,000 and 

□ 

sj3 n=512). Note that, unlike the case of FIG. 12, the smallest eigenvalues are much smaller than the 

■y largest eigenvalues. This suggests that the context vectors for this case could be adequately 

.sss. 

i=5 represented in a lower-dimensional space. 

To see the effect of changing n on the eigenvalues of T we generated two sets of context 
20 vectors: one set for n=280 and one set with n=5 12. Once the context vectors were built (using a 
1,000 document test corpus of Wall Street Journal articles from the TIPSTER collection) we then 
formed the matrix T, calculated its eigenvalues, and plotted them from largest to smallest. The 
results for n=280 and n=512 are shown in FIGS. 12 and 13, respectively. Clearly, for n=512 we 
see that the space can be reduced somewhat in dimensionality using the S VD method. Thus, we 
25 anticipate that the SVD method will be of use in developing sets of context vectors with the 

smallest possible dimensionality, without reducing the dimensionality too much (as in the n=280 
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case of FIG. 12. For most text data access applications we believe that context vectors with 
dimensionalities between 400 and 1000 will be adequate. 

We now consider the value of the SVD method for applications of context vectors. 
6. Practical Context Vector Generation Using SVD 
5 The main incentive for developing an SVD method for context vector generation is that it 

will enable us to rapidly build highly accurate context vectors for a core set of element classes. 
These context vectors can then be used with a two-pass bootstrapping method to build the 
context vectors for all remaining element classes. By core set it is meant that there often exists a 
set of high-frequency-of-appearance elements that can be used as the foundation for the context 
20 vectors for the remainder of the classes. For example, in text we might choose, say 2,000 high- 
j «3 frequency stems to use as a core class set. We would then use a very large corpus to compute the 
^ mutual co-occurrence matrix for this core class set (R would be a 2,000x2,000 matrix). The SVD 

!"U method would then be used to create a reduced-dimensionality set of context vectors from this R 

! y 

LU matrix. 

ie 15 The benefit of this method over the current initial random context vector/bootstrapping 

; jj method is that only one pass through the data set would be required to build an excellent set of 

ru 

y core context vectors. Further, since the co-occurrence matrix need only concentrate on the 
J«3 specific classes found in the core set, the number of calculations required during the one pass 
through the data is greatly reduced (as opposed to the current method, where all of the class 
20 context vectors must be adjusted during each of multiple passes through the data). This will 
allow the SVD method to use a much larger data set in the construction of the R matrix-thus 
yielding more accurate context vectors. 

Once the R matrix is formed, it will be necessary to apply the SVD method to it. Since R 
will be a large, sparse, symmetric, positive semi-definite matrix we will apply special numerical 
25 techniques for this process. A number of these have been developed. For example, Professor 
Martin Schultz of Yale University has developed a large library of software subroutines for this 
purpose (they are sold by Scientific Computing Associates, Inc.). His software has been used to 
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calculate the eigenvalues and eigenvectors for matrices of the above kind with dimensions as 
high as 100,000X100,000. 

Once the eigenvalues of R have been calculated we will then use a set of rules to 
determine where to set the dimensionalities of the context vectors. The context vectors for the 
5 core classes will then be created. 

Another possible approach to the SVD method is to use an adaptive neural network 
method for developing the eigenvalues and eigenvectors of the R matrix. This method has shown 
great promise, as it allows calculation of these quantities with typically greatly reduced 
computational burden, when compared with the algorithmic methods discussed above (see [1,2, 
W 3, 7, 10,1 1] for details on these neural network SVD methods). The primary reason for this 
*Z reduction is that the R matrix does not have to be computed. These methods work by passing 

| through the data itself. However, these methods will have to be adapted to this application, as our 
J* requirements are quite different from those of other applications that have used these methods 

^ (e.g., data compression and grammar analysis). 

s 15 At least two different approaches to building a complete set of context vectors can be 

-—3 

J3 tried using the above core context vector generation procedure. One approach is to simply let the 

y 

core set include almost all the classes. For example, in a typical corpus of approximately 1 
5 million newspaper articles there would be approximately 40,000 stems which would appear more 
than 10 times in the corpus. If we restricted this to stems which appear 100 or more times this 

20 number would be reduced to roughly 10,000. It may be feasible to directly compute the R matrix 
for these classes and carry out SVD on it. The remaining stems (about 1 million, approximately 
75% of which appear only once each) can then have their context vectors computed by adding up 
the weighted core context vectors in their co-occurrence sets during a single pass through the 
corpus. In this case, we believe that the SVD method will provide an excellent set of context 

25 vectors that have their mutual geometric relationships arranged exactly as dictated by the data 
itself. 
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Another approach to the development of a complete set of context vectors is to use a 
much smaller core set (e.g., the 2,000 highest frequency stems in the case of the TIPSTER 
corpus). Once the core context vectors have been created, they will be frozen. These core context 
vectors will then be used to determine the other class context vectors. Two passes of the batch 
5 bootstrapping method will be used. The first bootstrapping pass through the training data will be 
used to collect weighted sums of core context vectors for each co-occurrence set which contains 
a sufficient number of core context vectors. Those non-core classes that have sufficient numbers 
of co-occurrence sets will have their context vectors fixed at the end of this pass. A subsequent 
second pass will most likely be sufficient to calculate the context vectors of the remainder of the 
10 classes. It may be desirable to then go back and complete one more pass to readjust the core and 
w first-pass context vectors based upon the second pass results. This method requires three passes 
£H through the data, but lowers the size requirement for the core class set. 
j 8 y Both of these S VD context vector generation methods are expected to be much more 

ru 

Lu efficient and accurate than the current random initialization plus bootstrapping method. 

j; 25 Significant improvements in retrieval and routing performance are expected from this new 

O 

. ? s method of context vector generation. 

! j; 7. Context Vectors for Foreign Language Text 

One of the pressing needs in data access for text is to be able to retrieve and route 
documents in all languages that discuss a particular topic described in a single language query. 
20 HNC has devised a specific method for solving this problem using context vectors. This method 
is described in the this section using the example of a corpus containing documents in both 
English and Spanish. 

The first step is to build (or adopt from another system) a set of context vectors for the 
English portion of the corpus. A person who is skilled in both languages in then engaged. Their 
25 job is to create a list of tie words, which are words that have the exact same meaning in both 
languages. In informal discussions with persons who know both English and Spanish, English 
and Russian, English and Chinese, and English and Japanese it is clear that it is easy to produce 
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tie word lists having hundreds of entries for each of these language pairs. With effort, lists of at 
least two thousand tie words could probably be produced. We presume that this will be possible 
for any two human languages. To make this process simple, all tie words would be between 
English and the foreign language (in the case of this section, Spanish). This would seem 
5 reasonable because English is arguably the richest human language and it is also second in 
worldwide popularity (behind Chinese). 

Once the tie words have been selected, their context vectors in English are looked up. 
These are then transferred to the context vector set for the corresponding words in the foreign 
language. These tie word context vectors are then frozen. The remainder of the context vectors 
w for the language (a stop list and a stemmer can be used, if desired, but we do not believe that 
'"■T these are all that beneficial) are then created, using, for example, a neural network SVD method 
\f\ that adaptively builds new context vectors using the frozen tie word context vectors as a 

substrate. The final result is a set of context vectors for the foreign language that lie in the same 
W context space as the English context vectors. 

« 25 Key Point: Context Space is a universal meaning representation domain that we expect 

:J3 will be usable for representing the meaning of essentially all human data items (text in all 

ru 

! y languages, speech in all languages, imagery, and video). 

□ 

q Once the context vectors for the new foreign language have been built then the 

documents in that language have context vectors built for them and these are logged into the 
20 context space database. 

In essence, context space becomes a crude universal method of describing the usage or 
meaning of a data item. The same space is used for all languages (and, eventually, for video, 
imagery, and sound as well). By means of quires expressed as vectors in this universal context 
space we can retrieve data in all languages and media on the basis of its content. When finally 
25 built, this will be the ultimate content addressable memory system. 

To retrieve or route documents in multiple languages requires no new mechanisms. A 
query in one of the available languages is first formulated (the machine must be told which 
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language it is). This query is then converted to a context vector. The context vector is then used 
to search through the document context vector database to find the closest matches, which are 
then presented in rank order (Boolean queries can also be used, but the key words will only be 
used with documents in the language from which they cdme-which must be specified). The net 
5 result is a list of the most relevant documents, independent of language. The user can then select 
which documents they care to look at and they will be displayed in the text window of the 
system. 

In the case of a English and Spanish system, several advantages combine to make such a 
system much easier to build than, say, an English and Japanese system. First, significant volumes 
10 of Spanish text are available on CD-ROM. Second, the characters used in Spanish already exist 
l :f within our Match Plus, system. Third, many expert speakers of both English and Spanish are 

til 

SCI 

" : readily available to us. For these reasons, we believe that it will be possible to build an English 
^ and Spanish system. 

y 

^ Another potential advantage of having a common context space for all languages is that it 

" — 3 

15 will probably also be possible to build a crude gisting system. The idea of this would be that a 
J3 foreign language document would be displayed in a text windows. The user would call up a 

t: ; 

**j gisting feature (by using a keyboard control sequence or by means of a mouse-activated menu 
3 selection). The gisting feature would place a window approximately one paragraph long in the 
text. Next to this window (to the side of the highlighted text in the window) would be a second 
20 window containing a selection of English words that have context vectors closely aligned with 
the aggregate context vector of the material in the foreign language window (which would be 
computed automatically by the gisting system). The English words in the gisting window would 
be presented in order of closeness to the context vector of the foreign language window. 
Although they would not be formed into sentences, we believe that these words would 
25 nonetheless give a very clear idea of the content of the selected passage. The user could then 
scroll the gisting window up and down the foreign language text to survey its content. Further, 
the existing Match Plus highlighting system could be used to locate those passages of the text 
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that are most highly related to the subject mane of the query currently presented in the query 
window. In the end, we believe that this gisting window feature will, in many cases, obviate the 
need for translation of what are later recognized as irrelevant documents into English. This is yet 
another analyst productivity enhancement that we expect to flow from context vector technology. 
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DOCUVERSE 

An Intuitive Visual Representation System for Exploring Large Document Sets 
1 . Executive Summary 

A critical problem faced by analysts is the ever growing volume of written material that is 
available. It is said that well in excess of 90% of all desired intelligence information is sitting in 
available documents, waiting to be found and digested. Boolean query based document retrieval 
and routing systems were, historically, the first attempt to find a way to access large document 
corpora on the basis of topics of interest. Next-generation meaning-similarity based document 
retrieval and routing systems (in particular, HNCs TIPSTER program system) are now being 
developed. These are expected to significantly increase the productivity of analysts in terms of 
increasing their ability to much more quickly and comprehensively access documents that pertain 
to a particular area of interest. 

Although high-performance document retrieval and routing is of critical importance, it 
only addresses one aspect of an analyst's job. Namely, to probe document databases for 
information of a known type. Another, equally important, analyst function is to search for 
unexpected and unusual developments. To look for new trends and emerging activity patterns. 
Document routing and retrieval systems based upon searching for specified types of information 
cannot, by their very nature, be of much use in carrying out this exploratory function. What is 
needed is a way for analysts to somehow "get into" the universe of documents and roam around— 
to get to know the structure of the document database and look for any oddities or surprises that 
present themselves. 

This White Paper proposes a new type of document database exploration tool—one that is 
designed to allow analysts to exploit their powerful natural visual pattern recognition skills to 
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explore the information content of huge numbers of documents simultaneously. It will be the 
analyses window into the document universe. We therefore call this concept the DOCUVERSE 
system. 

The goal of the DOCUVERSE system is to support an intuitive, easy to control, 

5 exploration process whereby aspects of the contents of large numbers of documents can be 

rapidly assessed. The substrate for this process is a color graphics visual representation of a set of 
documents on a computer screen (e.g., a color graphics workstation or a PC or PS/2 computer . 
equipped with a graphics board). This representation, which will exploit some of the artificial 
reality image generation technology developed for use in Hollywood films and flight simulators, 

20 will portray each document, in the set of documents being considered, as a 3-dimensional object 
or icon with a specific shape, size, color, texture, and movement located in a three-dimensional 
cyberspace. Each of these attributes associated with each document corresponds to the 
document's content's similarity of meaning with one of a set of user-chosen features (such as a 
body of text indicating a topic of interest). Coded information about a document (such as when it 

15 was written or the identity of its author) can also be used as a feature. This will allow an analyst 
to view large groups of documents in a multi-attribute space. 

The project discussed in this White Paper will develop a prototype DOCUVERSE 
system. The system will be written in the C software language and it will run under Unix on an 
off-the-shelf color graphics workstation. It will not depend upon any other special hardware, but 

20 will utilize the workstation's color graphics display capability and the workstation's mouse. 

This project will also include the procurement and installation of a TIPSTER Software 
Evaluation System Sun Microsystems workstation system (integrated with the DOCUVERSE 
system workstation) capable of running the software being developed by the TIPSTER document 
detection contractors. 

25 2. Definition of the Problem 

2.1 The End-User Problem 
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Effective access to large databases of textual information is a continuing operational 
problem. Ongoing developments in computer networks, query-based retrieval and routing 
systems, and electronic mail systems are providing ever increasing access to documents of 
interest on specific topics. However, analysis of the contents of these documents and exploratory 
5 discovery of trends and changes in this content must still be carried out manually. Tools for 
supporting this crucial area of work simply do not exist Since these activities represent perhaps 
half of an analyst's job, the development of computer-based tools in this area is of the highest 
importance. 

2.2 The Technical Problem 
20 Advanced searching and retrieval methods provide the capability to locate a large fraction 

!; :f of the information on a specified topic that exists within a document database. The final product 
j-H of a system based upon these methods is a prioritized list of documents that are relevant to a 
!"U specified topic. The analyst must then examine these documents and use their content to reach 

ry 

\aI conclusions. Thus, a "one-document-at-a-time" analysis bottleneck is created, which often limits 

a 15 the analyst's ability to quickly identify trends, changes, etc. What is needed is another way to deal 

□ 

.fir with the contents of a preselected set of documents (e.g., the output of a query-based retrieval 
i % j and routing system) on a mass basis. Specifically, it is desirable to have a capability for viewing 
!;~ the documents as individual objects floating in a visual cyberspace, with the position and display 

of each object determined by its ranking or correlation with user-defined document attributes. 
20 The goal is to produce a document set content exploration system that can operate on large 
document sets and which can exploit natural human visual scene analysis capabilities. 
3. The DOCUVERSE System Concept 

Human vision can perceive and interpret many dimensions of information if the 
information is encoded and represented correctly. Context vector representation and high 
25 resolution displays are enabling technologies for visualization of textual information. Display of 
information can be accomplished such that humans can "navigate" through abstract 
representations of textual databases. That is the purpose of the DOCUVERSE system. 

- 47 - 5029677 



This section begins with a review of the HNC technique of context vectors. The 
DOCUVERSE System concept is then explained via discussions of its constitutive elements. 
3 . 1 Context Vectors 

HNCs TIPSTER project document retrieval and routing system and the DOCUVERSE 
system proposed here are both based upon the use of context vectors. A context vector is a point 
on the unit radius sphere in n-dimensional Euclidean space that represents the meaning content of 
a document (or a segment of a document). The components of the vector are the correlations of 
the overall document meaning with the meanings of some fixed features. Features include 
carefully chosen highly descriptive terms that represent the concepts of a small context-free 
language. Other future values are adaptively learned from the corpus. By means of a 
mathematical technique (see the HNC TIPSTER Project Proposal and Project Documents for 
details) a context vector can be assigned to any word, phrase, or document segment. As we have 
shown on the TIPSTER project, these context vectors provide an accurate guide to similarity of 
meaning— regardless of the nature of that similarity. This new approach can be contrasted with 
more conventional similarity of meaning techniques, such as the WordNet system. 

Prof. George Miller and his colleagues at Princeton University have developed a type of 
associational dictionary called WordNet. WordNet consists of several data structures— one for 
verbs, one for nouns, one for adjectives, etc. The verb structure is a mathematical ordered tree 
structure generated by the relation "in the manner of. A verb in the tree is connected above 
another verb if the latter verb action is in the manner of the former. For example, the verb walk 
lies below, and is connected to, the verb move in the verb tree, because walking is an action in 
( the manner of moving. The noun tree in WordNet is similar to the verb tree, except that the 
relation is "is a kind of. Thus, animal lies above and is connected to dog, because dogs are a 
kind of animal. 

Experiments carried out under the TIPSTER Program have shown that HNCs context 
vectors capture the very limited kind of similarity found in WordNet, and much more. Although 
exhaustive testing has not been carried out, spot checking has shown that pairs of words 
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WordNet considers close in meaning are also considered close in meaning by the HNC context 
vector system. However, context vectors go beyond this. For example, because of the very 
limited relations that are coded in WordNet, no relationship at all would be found for the words 
drive and automobile (since one is a verb and the other a noun). However, the context vectors for 
5 these words are quite close, as are those for block and city, and block and tackle. On the other 
hand, the context vectors for the words automobile and dog are not close. In fact, since the 
similarity of context vectors for words flows from the statistics of their close proximity in huge 
volumes of text, essentially all types of similarity of meaning become automatically incorporated 
into this meaning representation. This carries over into the context vectors for phrases and 
10 document segments, since these are constructed by taking weighted vector sums of word context 
;~ vectors. 

;; ? ! Besides providing an excellent vehicle for the encoding of meaning, the mathematical 

! a H form of the context vectors can be exploited to develop a fast searching method so that the 

ru 

similarity of meaning of each document in a document set can be determined rapidly. It is also 
15 possible to do trimmed searches to find only those documents that have a high correlation of 

! fl meaning with a particular context vector. These properties flow from the facts that comparing 

i y 

1 j similarity of meaning is carried out by the simple mathematical operation of taking the inner or 

Q 

□ dot product between the selected context vector and those of the documents to be rated as to 

similarity. Searches through large document sets for close matches can be carried out by means 
20 of hierarchical cluster searches, which are very efficient and fast. 

In summary, HNC's context vector approach provides a practical means for representing 
the meaning of a word, phrase, or document and for comparing the similarity of meaning of 
multiple documents. 

3.2 Data Objects 

25 The documents that make up the databases of interest for analysis contain information 

that is hierarchically organized into sets of nested data objects (see FIG. 14) For the purposes of 
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this proposal, these data objects could be any of the following (other possibilities, such as 
images, line drawings, and logos will be ignored here): 
Words. 

Part of a document (sentence, paragraph, section or chapter). 
5 Document. 

Set of documents. 
Database of documents. 
Set of databases. 

The DOCUVERSE system will be able to display content information on most or all of 
20 these types of objects. However, to keep the presentation simple, we will concentrate here only 
; '3 on the case where the data objects being portrayed are all documents. 
\f I 3.3 Artificial Reality Display Systems 

[U The DOCUVERSE system will employ state-of-the-art computer graphics techniques to 

W create a visual display of the document universe to be explored. This display will be constructed 
.'s 15 using artificial reality software that will be adapted for this project from an existing software 
!..□ package. The basic ideas behind this software are described in this subsection. 

y During the 1980s, two very similar standard high-level graphics software interfaces were 

Q 

i»3 developed~GKS (the Graphical Kernel System) and PHIGS (the Programmer's Hierarchical 
Interactive Graphics System). These are now merged into a standard PHIGS interface that is 

20 supported by multiple workstation vendors (e.g., Sun Microsystems, Hewlett-Packard, IBM, and 
Silicon Graphics), along with many enhancements (which we shall automatically include when 
we hereinafter say "PHIGS"). For details on computer graphics, GKS, PHIGS, X- Windows, and 
related issues see: Hill, F. S., Computer Graphics, Macmillan, New York, 1990; Kessener, L. R. 
A., Data Structures for Raster Graphics, Springer-Verlag, Berlin, 1985; and Foley, J. D., and van 

25 Dam, Fundamentals of Interactive Computer Graphics, Addison- Wesley, Reading, Mass., 1983. 
The basic idea of PHIGS is that the graphics creation job is broken into two pieces: 
describing the objects to be displayed to the user, and generating the display of those objects 
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from a user-selected eyeball position, direction of look, and panoramic angle of view (i.e., zoom 
level). The first job is that taken on by the application software. The second job is automatically 
carried out by PHIGS. 

The description of a 3-dimensional scene used as input to PHIGS is simply the definition 
5 of each individual object to be displayed, expressed as a set of linked polygons located in a fixed 
3-dimensional coordinate system, with each polygon having specified light reflectivity properties 
(color, specularity, texture, etc.). The polygons make up the visible exterior surfaces of the 
objects to be displayed. The lighting of the objects and the calculation of their appearance to the 
user from his or her selected eyeball vantage point is the job of PHIGS. 
10 On some workstations, PHIGS is implemented primarily in software. On others, PHIGS 

le 2 is implemented primarily in hardware. While the ultimate screen appearance is essentially the 
;^ same in both cases, those workstations that implement PHIGS in hardware are usually able to 
j^j display imagery at speeds that would be consistent with a human analyst's work pace, whereas 
W those that implement PHIGS primarily in software would be too slow for such applications. It is 

» 15 important to note that PC and PS/2 computers could also run PHIGS at a satisfactory speed, if a 

Q 

i,y special graphics display board were used. 

ru 

\f The artificial reality software that we will use on this project will be able to take the icons 

;5 to be displayed (with all of the "attribute" information about them-see below) and convert their 
attributes into a list of polygons that PHIGS should display to the user. This sounds simple, but 
20 will actually involve considerable complication. For example, if a group of icons exists at a 

distance that is too far away to flow the individual icons to be resolved, the software will have to 
replace the individual icons with a realistic rendering of a "cloud" of icons. In order to allow the 
effective use of human visual scene analysis, this rendering, and all others, will have to be done 
with very high quality. Another example is how icons of different size will be displayed at 
25 different distances. This will require some clever darkening or shadowing to make size obvious 
at variable distances. Another possibility is to employ a "fog" that reduces visibility at greater 
distances. These are all problems that have been solved by artificial reality companies that do 
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computer-generated graphics for Hollywood films, flight simulators, molecular modeling, high- 
quality video arcade games, and advertisements. 
3.4 Data Object Representation 

The DOCUVERSE system will display data objects to the user in the form of three- 
dimensional graphics icons. These icons will have attributes associated with them. An attribute is 
a user-selected descriptive feature that describes one aspect of the meaning of the data object. 
Each attribute will have a value between some minimum (say, 0) and some maximum (say, 1). 
The goal of the display is to show the user all of the attributes of each data object for a large set 
of data objects in one grand visual display. To do this, we shall exploit almost every human 
visual (and computer graphics!) skill. 

The basic idea is that the icon that represents each individual data object will have many 
visual degrees of freedom (see FIG. 15). 

These degrees of freedom will be "orthogonal", in the sense that they can be set 
independently of each other. Some examples are: position in space, size, shape, color, distortion, 
surface texture, and motion. Perhaps, via audio imaging, sound could also be used. 

Since the computers used will have graphics capabilities, the presentation of each icon 
will be visually very sophisticated. Each individually resolvable icon will be portrayed as an 
easily identified object in its proper position and possessing its assigned attributes. The display 
will resemble more a photograph of a real scene than an artificial reality display. Except that the 
objects will be icons. By means of simple mouse commands the user will be able to navigate 
around in the 3-dimensional cyberspace. A user-selected window will be available to show the 
entire cyberspace (a cube) and show the user's present position and their current direction of look 
through the screen of the display. The user will be able to effortlessly change their position to any 
new (x,y,z) position and rotate their direction of gaze to any new azimuth and elevation (the 
viewscreen will always remain horizontal, i.e., parallel to the x-y plane, in its boresight rotation). 
While quick changes in position and angle will be possible, smooth motion will not-as this 
would require an expensive specialized processor. 
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Icons that lie at too great a distance from the viewscreen to be resolved will be 
represented as points of light, or, if they are in groups, as shaded regions— sort of like the milky 
way. This will allow users to assess the locations of distant data objects. Options such as the 
introduction of a "fog" that limits distant vision will also be available. 
5 At low zoom levels the user will be able to see the entire cubical cyberspace within the 

viewscreen. At these zoom levels individual objects will typically be invisible and only large 
groups of documents will be seen. The density of these groups will be represented via 3D 
shading. By moving the viewscreen around the cyberspace cube the 3 -dimensional layout of the 
entire document set will be visible. This will then suggest strategies for moving into specific 
w regions for a closer look. 

Another user option will be to simply double-click on any icon to open that document for 
^ inspection. This will cause a new overlay window to appear containing the text of the object and 
i'U presenting the TIPSTER-like highlighter profiles (see TIPSTER Project documentation for 

ru 

W derails) for the object (for each of the icon display attributes the user has chosen). By clicking 

'•i. i 

ii 15 anywhere on any of the highlighter profiles the text block displayed will instantly jump to the 
i.,3 designated passage. The ability to roam a document universe at warp speed and then instantly 

ru 

\l locate and jump to passages of specific content in a document can reasonably be expected to 
l»5 increase analyst productivity many times over. The difference between manual searches through 
sets of retrieved documents and exploring the same set with the DOCUVERSE system will be 
20 like the difference between the Wright Flyer and the Starship Enterprise. We hope to give each 
analyst the power to continuously monitor all available textual information in their area of 
responsibility, and still have time to enjoy life. 
3.5 Attributes 

The attributes that determine the degrees of display freedom of the icons are chosen by 
25 the user by making selections via a user interface on the workstation. Preliminary concepts about 
how these attributes will be selected and used are now discussed. Task 1 of this project will be 
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devoted to carefully defining the various methods that will be available to the user for selecting 
and using the attributes. 

Coordinate System Selection: Coordinate Words 

Spatial coordinates of displayed icons are specified by three context vectors (see FIG. 4). 
Context vectors can be from words ("coordinate words"), groups of words, documents or sets of 
documents. Coordinate words define a (non-orthogonal) basis space for viewing the projected 
information content of data objects in 3D. Example coordinate words might be "Terrorist", 
"Drugs", and "Money Laundering". Positions for display of icons will be computed based on the 
projection (dot product) of the data objects' context vector with the context vector for the 
coordinate word. It will also be possible to use other data object attributes as "coordinate words". 
For example, one coordinate axis could be used to display the date of production of each data 
object. Another could display whether the author was Fidel Castro or not, etc. (See FIG. 16). 

Information Content Display Control: Icon Words 

In addition to the coordinate words, the user can specify a set of "icon words". Icon words 
determine the additional information that will be displayed (i.e., as the shape, color, texture, size, 
and motion of the icon, see FIG. 17). Icon attributes will be determined by computing the 
projection of the context vector for the data object attributes with context vector for each icon 
word. One icon attribute will be associated with each icon word. 

Qualification of Display Objects 

Qualification of data objects will be allowed to help reduce visual clutter and information 
overload. Boolean qualification will be sets of key words to match. Icon attributes will be 
displayed only for those objects that are a hit. Probably the most common qualification approach 
will be to simply use a retrieval and routing system to select a body of suitable documents. 

Exhaustive enumeration of data objects will be allowed: e.g. Document 4, Document 
1544, Document 3, Database "set.info.lib.text" either as keyboard command or "point and click" 
with mouse. 
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Context vectors can be used to perform data object qualification via dot products and 
thresholding (similar to routing approach). Clustering can be used to associate data objects with 
similar information content. Only objects within a cluster will be displayed. 

Traversing the Space: Changing Point of View 

Selection of coordinate words or object-associated data defines the display basis space. 
Selection of icon words determines which kinds of information are displayed. Input from a 
pointing device (mouse, etc.) determines the viewpoint of the user in the 3D coordinate space. 
Changes in this viewpoint (and in zoom level and angle of view) will be carried out using 
standard, off-the-shelf computer graphics and artificial reality methods. 

All data objects at the current level of display hierarchy are displayed using 3D, 
perspective and hidden line removal. User can modify his viewpoint, coordinate or icon 
specifiers at any time to gain full insight into the information contained in the database (see FIG. 
18). 

Scope Control: Preventing Information Overload 

Data object qualification provides the first method of overload reduction. Hierarchical 
organization of the data objects provides the second method of overload reduction. The user 
selects a data object for examination of subordinate objects using pointing device: "Zoom-in" to 
lower level of detail. When zooming in on an object, all higher level, non-selected objects are 
erased from the display. We will consider a mode whereby each document is nominally displayed 
as an icon and where the user will be allowed, if desired, to dive inside that icon and examine the 
document's chapters, sections and paragraphs as sub-icons. Diving inside a paragraph might 
cause the text of that paragraph to be automatically displayed in a screen window. 

At the lowest level in the hierarchy the text of the data object is displayed. The user may 
ascend and descend the hierarchy to any desired level. Higher level objects (such as sets of 
documents) are restored when ascending the hierarchy. 

Automated Determination of High Information Coordinates for use as Attributes 
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Besides allowing the user to select attributes, the DOCUVERSE system will be capable 
of automatically selecting coordinates for a document set. These coordinates will be selected on 
the basis that they have the highest possible information content. For example, the simplest such 
coordinates would be the principle components of the document set co variance matrix. These are 
5 obtained by considering the data object's context vectors as points in feature space. These points 
form a cloud, with one point for each object. Principal component analysis simply fits the best 
fitting ellipsoid to this cloud (in a mean squared error sense). The longest perpendicular principal 
axes of this ellipsoid are then the principal axes of the data set (in many real-world cases only a 
few of the ellipsoid axes will be very long— the short ones can then be ignored, with little 
20 resulting error). The attributes would then be these principal components. 

While principal component analysis usually works well as a first-order approximation, it 
;;■ ; often fails to yield accurate representations. A generalization of the principal component method 
^ has been developed by Cottrell, Monro, and Zipser and extended by Hecht-Nielsen. This method 
;«M uses a neural network to fit a general curvilinear coordinate system inside the data set, in which 
15 the coordinate axes remain approximately locally orthogonal, but curve to fit the actual form of 

! :D the geometrical "cloud" of data in the set. This method is expected to yield attributes that have 

i ?y= 

4 the highest possible information content in terms of their ability to represent the original 
q document context vectors. HNC is currently engaged in an effort to refine this method. The 
results of this effort, which is sponsored by the SDIO Office of Innovative Science and 
20 Technology and managed by the Mathematical Sciences Division of the Army Research Office, 
will be used in the project proposed here. 
3.6 A Specific Example 

Finally, we present a specific fictional example of how some of the capabilities of our 
proposed DOCUVERSE system might be used. Rather than present an example that explores the 
25 more obvious projected capabilities of the DOCUVERSE concept, we discuss a situation that 
illustrates one of the innovative ways in which the system will probably be used. The objective in 
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this example is to determine, in near real-time, the nature of the information being transmitted on 
a high-speed non-encoded (or decoded) data link. 

Data: Intercepted high bandwidth (10 MB/sec) data link that contains text information. 

Obstacles: There is far too much information for the available analysts to manually read 
5 messages within an acceptable amount of time. 

Approach: 

Step 1 : Specify Coordinates 

Analyst defines display coordinate space by selecting a set of coordinate words of 
interest: 
10 Drugs 

Money Laundering 
Terrorists 

DOCUVERSE displays representation of information contained in messages as 
icons in selected coordinate space (see FIG. 19). 
25 Step 2: Inspect and Select Data 

Analyst determines which messages are of interest in the "terrorist, drugs, money 
laundering" coordinate system by visual inspection and qualifies (selects) a set of messages with 
mouse for more detailed inspection (see FIG. 20). 
Step 3: Change Point of View 
20 Using a mouse, the analyst moves in the selected coordinate space to gain a clearer 

understanding of information in messages. The display is updated to reflect new point of view 
(see FIG. 21). 

Step 4: Select Icon Words and Attributes 

The analyst selects additional data dimensions for investigation by specifying icon words: 
25 "Sendero Luminoso" 

"Simon Bolivar" 
"Cinco de Mayo" 
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The display is updated to reflect the message information content relative to icon words. 
Messages that have high correspondence with these icon words will be easily visible (see FIG. 
22). 

Step 5: Zoom 

The analyst will then select the most relevant data object for zoom-in examination of 
message information. The display is updated to reflect adjusted point of view (see FIG. 23.) 1 1). 
Step 6: Inspect 

The text of the selected message is then displayed for manual inspection (see FIG. 24). 
3.7 Advantages of the DOCUVERSE System Concept 

In summary, the DOCUVERSE system will provide new data exploration capabilities for 
analysts. Some of the advantages of this system are listed below. 

High Speed: Allows inspection of the information content of large volumes of 
isolated text without the need for manual reading of messages. 

Cost Effective: Allows fewer analysts to perform more work. 

Easy to Use: Intuitive user interface exploits natural scene interpretation 
capabilities of human visual system. 

Compatible: DOCUVERSE uses same context vectors and features as HNC f s 
TIPSTER routing/retrieval system. 

Flexible: No assumptions made about nature of information. Could be applied to 
interpretation of speech if words are tagged with context vectors. Could even be used for 
specific-attribute gisting of foreign language message traffic. 

Standard Hardware: Standard computer graphics workstations or PC/PS2 
computers with a graphics board. 

Fast User Training: Intuitive human interface minimizes training requirement. 

sfc s(t ?|c sfc 



-58- 



5029677 



HNC MatchPlus Functional Specification 
HNC MatchPlus Modules 

The MatchPlus approach to information retrieval uses high dimensional vectors (280 
dimensions) called Context Vectors. Context vectors represent all stems, documents, and queries 
5 in the same vector space. This allows direct comparison between any two of these objects by 
taking a dot product. 

System generation begins with the preprocessing of documents to handle format problem, 
remove stopwords, etc. Next an inverted index is generated for later use in retrieval. The 
Bootstrap step creates context vectors for all stems, so that similarly used stem have similar 
10 vectors. Next, there is the option of generating word senses for commonly used stem, where each 
f «3 word sense receives its own context vector. The final step in system generation is the creation of 
;5 document context vectors by summing stem context vectors (with appropriate weightings). 
i«! For routing queries where there are documents with relevance judgements, a routing data 

; ^ generator extracts data so that neural network algorithm can generate a query. There are two 
~ 4 15 types of routing query, either term weights for a given query or an entire query context vector. 

i: 

^ For ad hoc queries, either the user supplies an initial query or the Topic Processor 

j y automatically converts a Tipster topic to an initial query. 

□ Once the query has been determined, the Document Retriever fetches as many documents 

□ 

as desired, in order by estimated relevance to the query. 
20 As a final step, the Output Combiner may optionally combine several lists of documents 

retrieved for the same query into one merged list. 

The following sections give detailed descriptions of MatchPlus modules in the same order 
in which they appear in FIGS. 25-28. 
1 Document Preprocessor 
25 1.1 Functional Description 

This module preprocesses raw text, taking stems off of words, taking out stop words, and 
deleting extraneous characters. 
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1.2 Data Specifications 

1.2.1 Input Specification 
Inputs consist of: 

The 'raw documents file' in the files of the corpus. 
The 'stopword file" which lists words to be dropped. 
The 'exception file" which lists words that are NOT to be stemmed. 
The 'word group file' which gives pairs or groups of words that are to be treated 
as a unit during training; an example is 'united states'. 
The "corpus descriptor file' which contains the location of corpus components. 

1.2.2 Output Specification 

Outputs consist of: 

The 'preprocessed documents file' which contains stemmed words, word pairs, 
etc. 

2 Inverted Index Generator 

2.1 Functional Description 

This module makes a pass over the preprocessed documents and creates the inverted 

index file and a file of information about each stem in the corpus. 
There are two steps in this process: 

1 . For each stem, the generator creates an inverted index consisting of the stem 

and the document(s) that stem occurred in. 

2. For each stem, the generator creates the 'stem info file" containing information 

such as the number of documents the stem occurred in, a pointer to the 
context vector for the stem, and the word-sense meaning of the stem. 

2.2 Data Specifications 
2.2.1 Input Specification 

Inputs consist of: 

The "preprocessed documents' containing stemmed words and word pairs. 

- 60 - 5029677 



2.2.2 Output Specification 
Outputs consist of: 

The 'inverted index file' consisting of stems and the documents they appear in. 
The "stem info file" consisting of information about stems such as number of 
5 documents the stem occurred in. 

3 Stem CV Generator (Bootstrap) 

3 . 1 Functional Description 

This module generates an initial random Context Vector for each stem and modifies each 
stem according to its neighboring stems. 
W Two passes are made through the entire corpus. For each stem in each document the three 

stems on either side (neighbors) are summed up. After the pass is completed the sum of all the 
neighbors for each particular stem is factored into that stem's context vector. 

3.2 Data Specifications 

3.2.1 Input Specification 
15 Inputs consist of: 

'Preprocessed documents' are needed for getting each stem and allocating a 
Context Vector (280 floats) and determining the stem targets and 
neighbors. 

The 'system op info file' specifies the document range in which bootstrapping is 
20 to occur. 

3.2.2 Output Specification 

Outputs consist of: 

The 'stem info file" with the trained stem CV's. 

4 Word Sense Generator 

25 4.1 Functional Description 

This module generates word sense information for frequently used stems. 
There are three main steps in this process: 
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1. For each appearance of a frequent stem, a window context vector is computed for 

surrounding stems. For space considerations, the number of windows for each 
stem is limited to about 300 randomly selected windows. 

2. The windows for each stem are clustered. Initially a variant of k-means clustering is 

used where the number of clusters is fixed at 7. Each cluster is represented by a 
centroid CV. 

3. For each cluster centroid, the closest 10 stems are found. This information is only 

useful for communicating with humans and for debugging; it is not used in 
retrieval. 
4.2 Data Specifications 

4.2.1 Input Specification 

Inputs consist of: 

The "stem info file' after bootstrapping. This file contains vectors for all stems. 
'Preprocessed documents' are needed for finding window CV's. 
The "system op info file' specifies which part of the text is to be used for 
collecting windows. 

4.2.2 Output Specification 

Outputs consist of: 

The "word sense file' specifying CV's for every word sense. 
The 'word sense stem info file', a modified stem info file that contains a pointer 
into the 'word sense file' for every stem that contains word senses. 
5 Document CV Generator 

5.1 Functional Description 

This module takes the preprocessed documents, stem information, and word sense 
information and creates document context vectors. The context vector for each stem in a 
document is weighted according to various formulae, summed together, then normalized to 
create the document CV. 
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5.2 Data Specifications 

5.2.1 Input Specification 

Inputs consist of: 

The 'preprocessed documents' containing stemmed words and word pairs. 
The 'system op info file' which contains information about the range of 

documents to process 
The "stem info file' which contains information about how many documents a 

stem occurs in, a pointer to its context vector, etc. 

5.2.2 Output Specification 
Outputs consist of: 

The 'document context vector file' which contains a vector for each document 
derived from the stem context vectors in that document. 
6 Routing Data Generator 

6.1 Functional Description 

This module generates the data necessary for stem routing and CV routing. 

1 . Stem routing data consists of the dot product between each stem in the query and the 

judged document. In addition the relevance judgement (a x 0' or T) is entered. 

2. CV routing data consists of the document's context vector and the relevance 

judgements. 

6.2 Data Specifications 
6.2.1 Input Specification 

Inputs consist of: 

The Query stems and their CV's. 

The 'relevance judgements' file with the official document id and the relevance 
judgements. 
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The 'read doc hash table' which is used to convert the official document id from 
the relevance file to an internal document number. This is necessary to get' 
the document CV. 

The 'context cvs doc file' which contains the document CV's. 
5 6.2.2 Output Specification 

Outputs consist of: 

The 'stem routing data file' which is an ascii file with the query stems and the dot 

products between the query stems and the judged documents. 
The 'CV routing data file', which is an ascii file with the products between each 
W query C V vector element and each judged document CV element. 

•J 7 Query Generator (Learning Algorithms) 
;f ! 7. 1 Functional Description 

This module generates a routing (or ad hoc) query based upon data generated by the 
^ Routing Data Generator. 

il 15 The learning algorithm used is the 'pocket algorithm with ratchet', a single-cell neural 

B network algorithm. 

j'fj 

l ---4 7.2 Data Specifications 

q 7.2.1 Input Specification 

For a Stem Routing Query, the inputs are a set of training examples, where each 
20 example consists of the relevance judgment for a document, plus the dot 

product of the document's context vector with the stem context vector for 
each stem in the original user query. 
For a CV Routing Query, the inputs are a set of training examples, where each 
example consists of the relevance judgment for a document, plus that 
25 document's context vector. 

7.2.2 Output Specification 

There are two types of query produced by this module: 
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A Stem Routing Query consists of a set of weights for each term in the original 
user query. 

A CV Routing Query is a context vector that can be used directly to take dot 

products with document context vectors in order to determine the order in 
which documents are retrieved. 

8 Topic Processor 

8.1 Functional Description 

This module generates the following from a Tipster Topic: 

1. The stems from which to calculate the query CV. The stems may be taken from 

any of the sections of the topic (e.g. concepts, definitions etc.). 

2. The stems that are required to be in a document (i.e. match terms). 

3. Overall weights to assign a particular section of the Topic. 

8.2 Data Specifications 

8.2.1 Input Specification 

Inputs consist of: 
An ascii Topic file. 

User inputs indicating the match criteria which sections of the Topic to use and 

any weights to apply to a particular section. 
The stop word list. 
The exception word list. 
The word group list. 
Word sense stem info file. 
Word sense file. 

8.2.2 Output Specification 

Outputs consist of: 
A user query. 

9 Query Processor 
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9.1 Functional Description 

This module generates a Query CV from a list of words taken from the following: 

1. Output from any one of the sections in the Tipster Topics (e.g. concepts section). 

2. Entered stems from the GUI 

9.2 Data Specifications 

9.2.1 Input Specification 

Inputs consist of: 

An ascii file with the query terms. 

The stop word list. 

The exception word list. 

The word group list. 

Word sense stem info file. 

Word sense file. 

9.2.2 Output Specification 

Outputs consist of: 
A Query CV (280 real numbers). 
10 Document Retriever 

10.1 Functional Description 

This module takes the output from the Query Processor and produces an ordered list of 
documents. 

First those documents that satisfy the SMatch criteria are identified. (This group may be 
empty if, for example, no document satisfies the criteria or there are no SMatch 
criteria.) These documents are then ordered by closeness of their context vectors 
with the Query Context Vector, and they are retrieved in that order. 

After the SMatch group is exhausted, remaining documents are ordered by closeness of 
their context vectors with the Query Context vector, and they are retrieved in that 
order. 
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10.2 Data Specifications 

10.2.1 Input Specification 

The Query Processor determines the SMatch criteria, as well as a Query Context 
Vector. 

5 The 'system op info file' specifies whether all or a subset of documents are 

eligible for retrieval. 
The 'inverted index' allows quick determination of those documents that satisfy 
the SMatch criteria. 

10.2.2 Output Specification 

10 Outputs are an ordered list of documents and dot product scores with the Query 

Context Vector. 
; 11 Output Combiner (Optional) 

11.1 Functional Description 

This module combines output lists, from MatchPlus or any other source, making use of an 

"■4. 

15 estimate of the quality of each source. 
=iD The algorithm is to give each document a merit score. Letting i run over all input lists, 

vj each document receives a score consisting of the sum of the quantity 'goodness of list i'/'position 
□ of document in list i'. 

Documents are then produced in order of their merit scores. 
20 11.2 Data Specifications 

11.2.1 Input Specification 

Several lists of documents. 

An estimated goodness for each list consisting of a non-negative fractional 
number. 

25 1 1 .2.2 Output Specification 

An ordered list of documents. 
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Section n. Image Retreival Using Context Vectors 

In this section n, all reference numbers refer to Figs. 29-43. 

FIG. 29 shows a block diagram of software and hardware components for implementing 
one embodiment of the present invention 100. Processor 102 is a conventional engineering 

5 workstation or other computer processor and may also be coupled to other processors accessible 
over conventional communications channels or buses (not shown). Processor 102 is 
conventionally coupled to storage 104 which may be a magnetic disk storage, a CD storage unit, 
or other conventional computer data storage unit. Storage 104 may also be coupled to other 
storage units accessible over conventional communications channels or buses (not shown). 

20 Storage 104 includes image database 106. Image database 106 contains images in electronic or 
computer-readable form. Images may be discreet images, video, or temporally related images 
(e.g., motion picture frames). 

Processor 102 is also conventionally coupled to memory 108 which is a random access 
memory (RAM) unit or other conventional computer memory. Items in memory 1 08 may 

15 alternatively be stored in storage 104 and accessed by processor 102 when required. Memory 
108 may comprise feature vector generation 1 10, prototypical feature vector generation 1 12, 
iterative training or context vector generation 1 14, indexing 116, clustering 1 17, and retrieval 
118 program components or modules as discussed below. Input 101 comprises conventional 
input devices such as a keyboard, mouse, trac ball, or touchscreen. A conventional display unit 

20 120 may also be conventionally coupled to processor 102. 

The preferred embodiment of the present invention may be implemented on any platform, 
operating system, and user interface of sufficient power and flexibility, such as: 
XWindows/MOTIF; Sun/OS SunView; Microsoft Windows, VAX/VMS, and the like, and may 
be operated in a distributed data / distributed processor environment, if desired. In the preferred 

25 embodiment, the invention is capable of handling a very large data corpus, e.g., containing over 
10 million images. 
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This aspect of the present invention is based upon an image representation scheme using 
context vectors. Context vectors are used to identify images of interest in existing image data 
bases and/or the routing and dissemination of images of interest in a real-time environment. A 
context vector is a fixed-length vector having a plurality of component values or weights that are 
determined based on relationships between meanings of information items or elements. In the 
present invention, context vectors represent the relative meaning among, or content of, 
information elements (described below) that are used to define images. The meaning is derived 
from geometric (e.g., proximity and co-occurrence of information elements in images) and 
quantitative statistics, data, and analysis (e.g., frequency of occurrence of information elements in 
the image database). Images are defined using information elements. Context vectors associated 
with each information element in an image are combined to form a summary vector. Similar 
images have closely aligned summary vectors, while dissimilar images have roughly orthogonal 
summary vectors. This representation scheme allows quantification of the similarity of images to 
be assessed by performing a simple dot product operation on the associated summary vectors, 
wherein a high dot product indicates a similar meaning. 

Note, the absolute orientation of a particular context or summary vector in the vector- 
space is irrelevant, as long as the relative orientation (with respect to other vectors) is 
representative of relative proximity of meaning and content. In other words, the problem of 
finding a set of vectors defining relative meaning has an infinite number of physical solutions in 
vector-space (absolute orientations) but only one logical solution (relative orientations). The 
context vector and summary vector generation scheme of the present invention produces a logical 
solution without regard to the absolute orientation of the vectors. 

Context vectors may be associated with words, terms, documents, document portions, 
queries, quantitative data, or any other type of information item or element. Many techniques for 
text-based context vector generation are applicable to image processing as well. See, e.g., 
Section I, Context Generation and Retrieval, above. 
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For image processing, context vectors are associated with certain image feature vectors. 
Elements in the feature vector, called features, are derived using affine wavelet transformations 
on sample points in images in the image database. Features are combined to form a feature 
vector. Such feature vectors provide an "image vocabulary" analogous to words in text. Due to 
5 the large number of possible feature vectors, a set of prototypical feature vectors, or atoms (or 
image information elements), are derived to form an "atomic vocabulary." Random high- 
dimensional context vectors are assigned to each atom. The context vectors are then modified 
according to the spatial relationship and co-occurrence of the atoms in the images in a procedure 
called bootstrapping. Compensation for frequency of occurrence of each atom in the database of 
10 images can also be made. Summary context vectors for each image are generated by combining 
the context vectors associated with the atoms that comprise the image. 

Context Vector Generation 

Referring now to FIG. 30, there is shown a high-level flowchart of the present invention. 
25 A corpus of images in image database 106 provides input of images 201 in electronic form to the 
system. The system is initialized 203 to generate atomic vocabulary 205 and stop list 207. 
Atomic bootstrapping 209 is then employed to produce context vectors 211 associated with each 
atom in the vocabulary. Context vectors 21 1 are used to generate 213 summary context vectors 
for each image 215. This method is described in more detail below. 

20 

Atomic Vocabulary Generation 

Referring now to FIG. 31, there is shown a flowchart of the process of generating atomic 
vocabulary 205 according to the present invention. Atomic vocabulary 205 is generated by using 
a combination of wavelet transformations, self-organization, and conventional statistics as 
25 described below. Image features, for example, primitive shapes, textures, or other similar 

attributes can be identified by wavelet transformations. These features are generally in a specific 
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spatial orientation to one another in the image. The features, their number, characteristics, and 
spatial orientations form the context and content information of the image. Thus, features in 
images are analogous to words in a document. Images can include any number of individual 
features in any combination. The information content of an image is defined in terms of specific 
features in a specific spatial orientation and position. When combined in the correct sequence 
and in the correct spatial positions, the features define the image. 

For each image 201 in image database 106, at least one wavelet transformation 301 is 
made of the image. Wavelets and wavelet transformations 301 are commonly known and are 
used to provide an efficient approach to localized characterization of information contained 
within images (see, for example, O. Rioul and M. Vetterli, "Wavelets and Signal Processing", 
IEEE SP MAGAZINE, October 1991; and, F. Hlawatsch and G.F. Boudreaux-Bartels, "Linear 
and Quadratic Time-Frequency Signal Representations", IEEE SP MAGAZINE, April 1992). In 
particular, the Gabor wavelet transformation provides a robust representation scheme of features 
that are sensitive to orientation and localized spatial frequency content. Thus, Gabor wavelets 
provide a mechanism for characterizing the information content of images in a compact and 
efficient fashion. 

In the present invention, multiple wavelet transformations centered at "sample points" on 
the image 201 are made 301. To implement this process, Gabor wavelet parameters are selected 
to capture orientation-sensitive spatial frequencies at a variety of orientation/ frequency 
combinations over the range of interest. Both sine and cosine wavelets are used to provide phase 
discrimination. For example, a combination of eight orientations, five spatial frequencies, and 
both sine and cosine (even and odd order) phases can be used. Thus, a total of 8*5*2, or 80 
wavelets are associated with each sample point on the image 201. The actual number of 
combinations for a particular application is determined empirically based on the detail and 
accuracy of image representation required and processing capabilities. Due to the high level of 
processing power required during wavelet transformation, the preferred embodiment can employ, 
for example, the Vision Processor (ViP) chip set to perform high-speed image processing 
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operations. The ViP chip set is available commercially from HNC, Inc., San Diego, CA, as a 
plug-in circuit board for workstations manufactured by Sun Microsystems, for example. 

A simple uniform grid-point approach can be used to determine the "sample points" on 
the image where the wavelet transformations are centered. Alternatively, more elaborate 
conventional "fixation control" techniques to determine "points of interest" within images can be 
used (see, for example, Abbot, Lynn, "A Survey of Selective Fixation Control for Machine 
Vision" IEEE Control Systems, August 1992). The accuracy of the image representation is 
linked to both the number of points sampled and the discrimination of the wavelet 
transformations employed (i.e., the number of coefficients computed at each sample point). In 
the example above, 80 real values are associated with each sample point. These values comprise 
a "feature vector" 305. Atoms in an atomic vocabulary 205 are formed from a subset of the 
feature vectors 305 in the manner described below. 

The result of the wavelet transformations 301 and feature vector generation for each 
sample point on each image can result in an extremely large number of unique feature vectors 
305. For example, using only a binary representation for the wavelet transformation values could 
result in 2^0 (10^4) unique atoms. Accordingly, the present invention recognizes that a small 
fraction of the feature vectors 305 can sufficiently represent the images in the image corpus. 
Thus, the set of generated features vectors can be reduced to a smaller number of prototypical 
feature vectors using statistical and neural network techniques. 

In one embodiment, a vocabulary size, N, is chosen which defines the problem size 307. 
A reasonable range for N has been found to be between 10^ and 10^. The exact vocabulary size 
is empirically determined based on processing power available, the level of detail required for 
image definition, and the size and diversity of image database 106, for example. Statistical 
prototype feature vectors are selected using neural network self-organization techniques 309. 
Basically, the best N feature vectors 305 that represent the entire population of vectors are chosen 
by using a vector quantization (VQ) training operation. In the present invention, the VQ "learns" 
statistically representative exemplars given a population of feature vectors 305. Conventional 
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VQ training techniques 309 such as extensions to clustering techniques (such as the widely used 
K-means algorithm) or variations on self-organizing approaches ("Self Organization and 
Associate Memory" by Tuevo Kohonen, 2nd Ed, Springer Verlag, 1988 P. 199-209) can be 
used to transform a potentially infinite number of examples into a finite (but possibly large) 
5 atomic vocabulary of statistically representative examples 205. See, e.g., Section I, Context 
Generation and Retrieval, above, which describes a singular value decomposition (SVD) method 
for reducing the vector space. 

The prototypical feature vectors form atoms or information elements in the atomic 
vocabulary 205. Each atomic vocabulary item 205 is given a unique atom ID 31 1. As described 
10 below, each atom ID 31 1 so defined is then associated with a high-dimensional context vector to 

^ discern the proximal meaning between such atoms. The resulting VQ network is also saved 313. 

rl In one embodiment, further statistics of the entire set of feature vectors 305 are 

determined to generate a "stop list." As in textual applications, images have a significant number 
of high frequency-of-occurrence atoms that convey little, if any, information. In the textual 

^ 15 application analogy, words like "a", "an", "but", and "the", are put on a "stop list" and not used 

] S to defined the textual document. Similarly, statistics can be used to define high-occurrence, low 

i -sr 

%f informational atoms in images. 

□ Referring now to FIG. 32, there is shown a flowchart of the process of generating stop list 

207 according to the present invention. Feature vectors 305 are fed through the previously saved 

20 VQ network 313 and are thus quantified or mapped 401 to the nearest-neighbor equivalent 

atomic vocabulary vector 205. Atomic IDs 311 associated with the mapped atomic vocabulary 
vector 205 are thus generated 403. Frequency-of-occurrence statistics are computed and 
accumulated 405 for the atomic IDs 403. Atomic IDs 403 having a frequency exceeding a cutoff 
407 are placed on "stop list" 207. The cut-off frequency value 407 can be empirically 

25 determined based on the image database subject matter and level of image representation detail 
required, for example. Stop list 207 is used to eliminate low information atoms in when 
representing images using the atomic vocabulary 205, 311.. 
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Image Characterization and Representation 

Any image may be represented in terms of the unique set of information elements in the 
atomic vocabulary 205, 311. Referring now to FIG. 33, there is shown a flowchart of a method 
5 for representing a new image using atomic vocabulary 205. As discussed above in association 
with feature vector generation (see, 301, 305), feature vectors for a new image 501 are generated 
by performing a plurality of wavelet transformations at a plurality of sample points 503 to 
generate a feature vector 507 at each sample point in the new image 501 . As in the example 
above, each feature vector 507 may comprise a vector of 80 real values. 
10 The previously saved VQ network 3 13 is used to map 509 each feature vector 507 to the 

^ nearest-neighbor equivalent atomic vocabulary vector 205. This provides a good approximation 
!;*! of the original image 501 and constrains the number of possible combinations of feature vectors 

; B f; 507 used to represent the image 501. Atom IDs associated with each mapped atomic vocabulary 

iy 

J s H vector 205 are thus generated 511. Atoms IDs 511 contained in stop list 207 are removed as they 

"'"4 

; ;_ 15 add little information. Thus, an unlabeled, unconstrained image in electronic form 501 is 

transformed into an image representation 513 comprising atoms 511 from atomic vocabulary 

S! 205,311. 
□ 

□ The atomic representation 513 of image 501 provides a significant data compression over 

conventional image representation schemes. In fact, since the atom vocabulary 205 is relatively 
20 rich, the atomic representation 513 can often substitute for the original electronic image 

representation 501 for many purposes. An image can thus be represented by a set of codes 511, 
as defined by the atomic vocabulary 205, together with an indication of the location to each atom 
in the image (e.g., locations taken from the sampling grid used). 
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Context Vector Generation (Bootstrapping) 

The context and relationship of atoms 511 in an image and in the entire image database 
106 is learned from examples. This learning procedure is referred to as "bootstrapping," and is 
basically a type of self organization. The bootstrapping procedure is based on the concept of 
atom spatial and co-occurrence relationships. That is, atoms that appear in proximity to one 
another in a significant number of images are conceptually related (just as are words that appear 
in proximity to one another in a significant number of documents are conceptually related). The 
closer the atoms are, the stronger the relationship between them. Additionally, the relative 
importance of the atoms (and their influence on neighbors) is related to both the local frequency 
(within an image) and the global frequency (within the corpus of images). The specific technique 
for implementing these relationships is described below. 

Referring now to FIG. 34, a method for initializing the system before bootstrapping is 
shown. According to the present invention, each atom ED 3 1 1 (corresponding to atoms in atomic 
vocabulary 205) is associated 605 with a high-dimensional context vector 601. A high 
dimensional vector space is advantageous in that 1) initial values for each vector can be 
generated by using a random Gaussian distribution which results in nearly orthogonal vectors 
(indicating no initial relationship between context vectors); 2) the information capacity of the 
vector space is very large when real valued vectors are used; and, 3) vector superposition allows 
summation of independent vectors into a single vector and later recovery of independent vectors 
via dot product. Initial values of each element in the context vectors 601 are randomly generated 
by using a random Gaussian distribution 603. This results in nearly orthogonal vectors indicating 
no initial relationship between context vectors. Each atom ID 31 1 is randomly associated 605 
with an initialized context vector 603. 

Atomic bootstrapping 607 is then performed by a constrained optimization technique 
where an objective function contains terms for local atom relationships, as well as global terms to 
insure that the statistics of the entire vector set remains well-behaved. In the preferred 
embodiment, a commercially available SIMD Neurocomputer Array Processor (SNAP) is used as 
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a hardware accelerator due to the compute-intensive nature of the learning operation. The SNAP 
hardware accelerator is available commercially from HNC, Inc., San Diego, CA. 

Referring now to FIG. 35, there is shown a flowchart of one embodiment of the 
bootstrapping procedure 607. Alternative embodiments that are optimized for processing speed 
5 may also be employed. In general, the bootstrapping procedure allows all atoms 511 within an 
image 501 to influence each other to varying degrees depending on proximity, and co- 
occurrence. The intra-image influence between atoms 51 1 is distance-weighted according to a 
Gaussian distribution in that atoms that appear close together in an image are given higher 
weight. In the preferred embodiment, all atoms 511 influence all other atoms within a single 
20 image, although this is not necessary. Alternatively, the influence could be limited to atoms 
' ! 1 within a certain "radius" of the target atom. This approach is similar to a Parzan window 

;f ; approach that performs non-parametric estimation of probability density based on superposition 

"'J 

jjjj of Gaussian weighted samples. See, e.g., Parzan, E. "On Estimation of a Probability Density 
J'H Function," Annals of Mathematical Statistics 33, 1065-1076, September 1962 and Rosenblatt, M. 

" s 15 "Remarks On Some Nonparametric Estimates Of A Density Function," Annals of Mathematical 
'•■■3 Statistics, 27, 832-837, 1956. The vector co-occurrence of atoms 51 1 is learned in proportion to 



their mutual spatial co-occurrence frequencies across all images 201 in image database 106. 
The following procedure is performed during each learning pass. Multiple learning 

passes through the corpus or database of images are generally required. First 700, the system 
20 selects 701 an image 501. Typically, images are selected sequentially, starting with the first 

image and proceeding through the last image. However, the learning approach makes no 

assumptions concerning the order of presentation of the images. For the selected image 701, the 

atoms 511, and associated context vectors (CVs) 605, that comprise the image are retrieved 703. 

In one embodiment, an image data structure contains a pointer to the linked list of atoms 511 that 
25 represent the image 501, and their associated context vectors 605. 

An atom 51 1 of the selected image is targeted 705. Typically, atoms are targeted 

sequentially, but this is not necessary. The associated context vector for the target atom 
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(Target_CV(old)) is retrieved 707. A zero-valued temporary context vector for local summation 
(Temporary_CV) is also created 709. 

Each atom 51 1 in the image 501 that is not the target atom 705 is selected 71 1 in turn, 
and the associated context vector of the selected atom is retrieved 713. Alternatively, each atom 
5 in a defined "learning region" or "radius" around the target atom 705 can be used (instead of all 
atoms in the image), for example, all atoms within 2 sample points from the target atom 705. 
The physical Euclidean distance, D (a scalar), between the target atom 707 sample point and the 
selected atom 713 sample point is computed 715 (e.g., points and distances defined in sample 
point grid). This distance is pictorially shown in FIG. 36. The weighted function, G, of D is 
1 o computed 7 1 7 as : 

G = e'^ /a (Eq. 14) 

where a is the standard deviation of a Gaussian distribution used in the proximity 
weighting function. The context vector of the selected atom 713 is weighted 719 by G and 
inversely by the frequency of occurrence of the selected atom 51 1 in the image 501 (an indicator 
^ 15 of importance), fi Alternatively is the frequency of occurrence of atom 51 1 in the entire image 
i'3 database 106. The weighted context vector is added 721 to the temporary context vector 709. 

The steps 711 through 721 are repeated 723 for each atom in the image that is not the 
target atom 705. Once all atoms 711 have been selected 723 (other than target atom 705), the 
system computes 725 an error vector (Error_CV) for the target context vector 707, according to 
20 the equation: 

Error_CV = Temporary_CV - Target_CV(old). (Eq. 1 5) 

ErrorCV may also be frequency weighted in proportion to the frequency of occurrence 
of the targeted atom in the corpus. Error_CV is squared and added to a variable Mean_SQ_Error 
725 to keep track of the total mean squared error metric for the image corpus. 
25 The system then updates Target_CV 727 according to the equation: 

Target_CV(new) = Target_CV(old) + a *Error_CV (Eq. 16) 
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where a is a predefined proximity constraint or step size that limits the correcting effect 
ofError_CV. The variable a is empirically determined based on variances in the image 
database and a desired learning rate, for example. Additionally, constraints may be added to 
prevent the target atom from "collapsing in" on the selected atoms. See, e.g., Section I, Context 

5 Generation and Retrieval, above. If there are atoms 511 that have not been targeted remaining in 
the selected image 729, the system returns to step 705. If not, it checks 731 whether there are 
more images in the corpus. If so, it returns to step 701; if not, the learning pass ends 733. The 
learning pass may be repeated multiple times. Generally, the process in FIG. 35 is repeated until 
the rate of change of Mean_SQ_Error 725 is arbitrarily small (i.e., the system is not learning 

w anything new). The cutoff value of Mean_SQ_Error is determined empirically based on 
variances in the image database and accuracy required. Once the context vectors 727 are 
iteratively determined, they are saved in storage 104. 



•-4 
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The procedure of FIG. 35 is expressed mathematically by the equation: 

cv;: gel = cc + G{\cvZ - cv\) • cv, ( Eq . n) 



target 



,0 

pi : where CVi is the selected context vector and fi is the frequency of occurrence of atom / 



in the image (or alternatively the entire image database), and G() is Gaussian weighting e ° la . 



Summary Vector Generation 

20 Context vectors 727 associated with the atoms 511 that comprise an image 501 are 

combined to form summary vectors for images. Referring now to FIG. 37, there is shown a 
flowchart of the steps that are performed for each atom 51 1 in the image 501, in order to generate 
a summary vector for the image 501 . First, 900, for each atom 901 in the image 501 , a check is 
made to determine 902 if atom 901 is on stop list 207. If atom 901 is on stop list 207, the next 

25 atom is processed 901 ; if not, the associated context vector 727 is retrieved 903. Context vector 
727 is weighted 905 by the function l/ln(N) where N is the frequency of occurrence of atom 901 
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in image database 106, 201 (alternatively, N is the number of images 201 that contain atom 901). 
Weighted context vector 905 is accumulated and summed 907. After all atoms 901 for image 
501 have been processed 909, the resulting summed vector is normalized 91 1 to produce 
normalized summary vector for the image 913. Normalization is performed by dividing each 
component in the vector by the absolute magnitude of the context vector. The magnitude of the 
context vector is determined by taking the square root of the square of all of the components in 
the context vector. 

Of course, alternative training techniques may be used. More generalized context vector 
techniques are disclosed in Section I, Context Generation and Retrieval, above. 

Summary Vector Storage (Clustering) 

Once the summary vectors 913 have been determined, they are stored in storage 104. As 
described above in Section I, summary vectors 913 can be stored to further reduce searching time 
by creating cluster trees. Referring to Fig. 38, an initial parent node 1001 at the top of a tree 
indexed as level 0, node 1, contains all of the summary vectors 913 for all images 201 in the 
image database 106. A series of child nodes 1005, each branching from the initial parent node 
1001, are created at a next level of the cluster tree by employing clustering algorithm 1003. For 
each child node 1005, further branching can be made by treating the child as the parent in 
clustering algorithm 1003. This process can be repeated until a sufficient level of clustering 
detail is achieved 1007, the result being centroid consistent clusters 1009. 

Fig. 39, further details a flowchart of clustering algorithm 1003. Centroid-consistent 
clustering algorithm 1003 divides parent summary vectors 1001 among a series of children nodes 
1005. A group of clusters is centroid-consistent if every member of every cluster belongs to the 
cluster in the group with the closest centroid. A centroid is determined by taking, for each 
dimension, an average of the component values from all of the summary vectors within the 
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group. In a preferred embodiment, the convergent k-means clustering algorithm, a well known 
clustering algorithm, is used. Convergent k-means clustering can be performed as follows: 

1 .Begin with any initial partition that groups the vectors into k clusters 1101. For 

example, take the first k summary vectors as single element clusters and set the 
initial value of the centroid of each cluster to equal its member vector 1 103. 
Assign each of the remaining summary vectors to the cluster having the nearest 
centroid (computed by dot product with each cluster) and recompute the centroid 
for the cluster which gains a vector 1 105; 

2. Take each summary vector in sequence and compute its distance from the centroid of 

each of the k-clusters 1 107 (computed by dot product). If the vector is not 
currently in the cluster with the closest centroid, move the vector to that cluster 
and update the centroids of the clusters that gain or lose a summary vector 1 109. 
This step is performed for each summary vector 1111. 

3. Repeat step 2 until convergence is achieved, that is until a pass through all of the 

summary vectors causes no new assignments 1113. Thus, centroid consistent 

clusters at a new level are created 1115. 
Since convergence 1113 may be rather time-consuming to achieve, the clustering 
algorithm can be simplified by limiting the number of repetitions of step 2 (1 107, 1 109, 1 1 1 1). 
After a given number of repetitions, such as 99, the centroids can be frozen. Then, one or more 
passes can be made through all of the summary vectors, distributing the vectors to appropriate 
clusters, but without updating the centroids. The centroids approximated in this manner will no 
longer be exact centroids, but such approximate centroids may be sufficient for most applications 
of the present invention. It is not necessary that the centroids be precise according to the present 
invention; however, it is preferable that the clusters be centroid-consistent with the approximate 
centroids. "Centroids" as used hereinafter in this application shall mean approximate 
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centroids — in other words, a centroid sufficient to establish centroid-consistent clusters. Each 

node is identified by its centroid for use in the searching process. 

In forming a next level of clusters 1007, the nodes in the level above become parent 

nodes to as set of child nodes below. Only the summary vectors assigned to a parent node are 
5 used in the clustering algorithm to form a the child nodes which branch from that parent. This is 

repeated across the entire level of parent nodes and on subsequent levels so that fewer and fewer 

context vectors are assigned to the child nodes on each lower level. 

The nodes form a tree pattern in which each node branches from a node in the level 

above. Thus, each summary vector is assigned to a node on each level of the cluster tree. Each 
10 node has a centroid. The bottom-level node assignments for each summary vector and the 

centroids for each node are stored for use in the search and retrieval algorithms. On the bottom 

level of the tree, each node points to each summary vector assigned to it. The nodes on the 

bottom level may be referred to as "buckets." 

Once a cluster tree has been set up 1009, it is a simple matter to add a new record 
15 summary vector to the tree. The initial branches of the tree are examined to find the closest 

centroid. The summary vector is assigned to the node with the closest centroid. Then, the 

branches from that node are examined for the closest child node centroid, and the process is 

continued until a bucket is reached. The new record is then assigned to the bucket with the 

closest centroid of those buckets branching from the node on the previous level to which the 
20 summary vector was assigned. The centroids are not changed. This action preserves centroid- 

consistency of the clusters. If a bucket gets too large, the summary vectors on the bucket can be 

divided into subclusters on a subsequent level. 

Summary Vector Indexing 

25 Images can also be indexed to terms. Index terms can be textual words or codes, for 

example. More than one index term can be associated with an image. For example, an image of 
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a dog may be indexed to the textual words "dog", "bark", and "pet". Thus, an index to textual 
terms allows for searching of images using textual retrieval techniques. The underlying principle 
is that images that have similar content will be associated with similar index terms. Index terms 
can also be associated with atoms or sets of atoms. The basic indexing algorithm as it relates to 
5 text is also described above in Section I. 

Referring now to Fig. 40, a high level data flow diagram of image indexing is shown. 
First, a set of pre-indexed "training images" is constructed, i.e., associated index terms are 
determined for each image 1201 . The index term/image association is generally derived by 
humans according to experience. One or more index terms may be associated with each image. 
20 An image context system is built according to the invention as described above 1203, i.e., an 
atomic vocabulary 205 with stop list 207 is constructed, bootstrapping is performed, and 
!;'! summary vectors 913 are generated for each image. 

LH A new image 1204 is characterized 1205 using the existing atomic vocabulary 205 and 

stop list 207 as described above (i.e., wavelet transformation and mapping to closest atoms in 

» 25 atomic vocabulary), and a summary vector 1206 for the image is generated. Similar images, and 

□ 

v3 their associated index terms, in the existing training corpus are retrieved by calculating a dot 

ru 

y product between the new image summary vector 1206 and each summary vectors 913 associated 
q with the images in the training corpus 1207. Optionally, a threshold can be set to retain the best 

matches, e.g., the top ten similar images and associated index terms 1208. The associated index 
20 terms retrieved are weighted by both similarity 1209 and frequency of occurrence 121 1 and an 

ordered list is produced of the most relevant index terms 1213. The new image is then associated 

with the listed order of index terms 1215. A threshold 1217 can also be used to choose the top N 

index terms from the weighted list. 
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Retrieval 

Referring now to Fig. 41 , there is shown a flowchart of the retrieval method of the present 
invention. The system allows the user to formulate image queries using several approaches. For 
example: 

•Image Queries: Whole images can be used as queries. This is, in essence, asking the 
system to retrieve images that "look like" a particular image. For this query mode, a 
summary vector for the query image is generated (i.e., using the atomic vocabulary 
205, 311) and is used as the query context vector. 
•Atomic Queries: The user can use a mouse or other computer input device to select 

portions of a displayed image to be used as a query. The system determines the atoms 
contained in the designated portions of the displayed image and forms a query context 
vector using the atomic vocabulary 205, 311. The query context vector is computed 
from a weighted sum of context vectors of the atoms in the portion of the displayed 
image being used as the query. 
•Textual Searches: As described above textual index terms or codes can be associated 
with images or atoms. For example, the context vector for the word stem "aircraft" 
could be aligned with a visual representation of an aircraft. Thus, text-based queries 
could be used in retrieval of image data. 
•Boolean Matches: The user may specifically identify atomic vocabulary codes to be 
matched by the search using boolean logic (e.g., and, or, not operators). Boolean 
logic may also be combined with, for example, the image query technique. 
Regardless of the approach used, user queries 1301 are converted into a query context 
vector 1303. For example, images and image portions are represented by the atomic vocabulary 
vectors 205 and a weighted query vector is formed from the atomic vocabulary vectors (see Fig. 
33). Image relevance is assessed by computing the dot product 1305 of each image summary 
vector 913 with the query context vector 1303, and accumulating the results 1307. The images 
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are sorted by dot product 1309. Images associated with the summary vectors that produced high 
dot products are retrieved 1311. 

If desired, image summary vectors may be organized in cluster trees for quicker 
searching. By using the above-described cluster tree storage mechanism (Figs. 38 and 39), the 
searching task can be greatly accelerated. The query vector is used to identify the centroid that is 
closest to the query vector (by dot product computation). The search is performed using a depth 
first tree walk. A branch is followed down the tree, taking the node at each level having the 
centroid closest to the query vector. The search proceeds down the tree until a bottom level node 
without children (e.g., a bucket) is reached. Each of the summary vectors in the bucket is 
compared with the query vector (again by dot product computation) to identify the closest 
summary vector(s). 

The system may also employ relevance feedback, whereby the user specifies which of the 
retrieved images are most helpful. A new search may then be performed using the summary 
vector for the specified images. This technique reduces the time required for searches and 
improves system effectiveness. 

Referring to Fig. 42, highlighting of retrieved images may also be provided. Highlighting 
indicates image portions having high correspondence to the query vector. To perform 
highlighting, the retrieved image is partitioned into blocks 1401. A context vector is computed 
for each block 1403, and the dot product of the query context vector 1303 and each block context 
vector is computed 1405. The block 1401 with the highest dot product has the highest degree of 
correlation with the query 1303, and is indicated by some visual distinction such as color, inverse 
video, or some other form of identification 1407. The highlighting concept is pictorially shown 
in Fig. 43. 

Also, as described above, summary vectors of images can be visually displayed for 
navigation through the corpus of images to find images that are of interest. Essentially, the 
summary vectors for images and other information items are displayed in a pseudo-orthogonal 
display having axes corresponding to query terms (the query terms need not actually be 
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orthogonal to one another). In this application query terms may be atomic vocabulary items, for 
example. Image summary vectors may also be displayed with visual attributes representing 
relative similarity of meaning with query terms, e.g., color, shape, orientation, texture. 

Practical applications of the present image storage and retrieval method and system are 
numerous. For example, the system and method can be used in medical imaging to find 
myopethis, tumors, or bone breaks in medical images and X-rays. Signature recognition of, for 
example, checks, or other valid signatures can be performed using images of such signatures in 
image database 106. 

As noted above, image database 106 can comprise video or temporally related images. 
Thus, video data can be used in the present invention. For example, films can be input into the 
system frame by frame (or using lower sampling frequencies, e.g., every 5th frame). Thus, films 
or portions thereof can be compared. 

As above, the use of context vectors to represent data from several different data domains 
in a common universal meaning space is possible (see Development of Context Vector™ by 
Singular Value Decomposition, above). In the textual application of context vectors, different 
languages are linked together by word-pairs or "tie" words that have the same meaning in both 
languages (e.g., water and wasser in English and German, respectively). The context vector of 
the tie word in English is used as the context vector of the word in the foreign language. 
Training of the context vectors in foreign language proceeds normally, except for the fixed 
context vector associated with the tie words. Context vectors in the foreign language that are 
close (e.g., dot product) to the tie word context vector will be similar in meaning to the English 
tie word. Thus, English queries can be used to retrieve foreign works. The joinder of foreign 
language data domains represents a common "universal meaning space." 

Similarly, a common universal meaning space can be constructed for images (or video) 
and sound, for example. Certain "ties" can be fixed between images (e.g., a St. Bernard dog) and 
audio clips ("bark") (context vectors in the audio domain may be, for example, frequency and/or 
amplitude measurements). Training of the context vectors in each domain proceeds normally 
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except for the fixed ties. Thus, a picture of a Chihuahua dog, would be highly linked to the 
sound clip "bark" by virtue of being close in meaning (i.e., summary vector dot product) to the 
St. Bernard dog image. Thus, audio and video data domains can be linked in a common or 
universal meaning space. 

Similarly, text and images can be linked, for example magazines and pictures in 
magazines. Here, the text surrounding the image (e.g., in a caption) can be automatically "tied" 
or associated with the image. The strength of the tie association can also be varied according to 
spatial proximity of the text to the image, boldness, font, or other factors. See, Section I, Context 
Generation and Retrieval, above, for an application of context vectors to text. Information 
elements in the textual application may be, for example, words or word stems (root words). 

The above description is included to illustrate the operation to the preferred embodiments 
and is not meant to limit the scope of the invention. The scope of the invention is limited only by 
the following claims. From the above discussion, many variations will be apparent to one skilled 
in the art that are yet encompassed by the spirit and scope of the invention. 
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